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Chapter 1. BASIC ANALOG BLOCKS 


1-1 Common Uses of the Analog Computer 

Although analog computers may be used for the solution of problems which do not Involve 
differential equations, they are often considered as differential equation solvers. For this 
reason, the analog computer has been called the "differential analyzer." 

As we know, differential equations are important in the study of dynamic systems because 
in these systems the Important variables are changing with respect to each other, so that their 
behavior may be expressed mathematically by differential equations. In some cases the physical 
variables are changing with respect to real time, and we are asked to find the so called 
"transient behavior" of the system. In other cases the variables are not changing with respect 
to time, and we may want to study the "steady-state" solution of the system. Although the 
analog computer may be used for steady-state problems as well as for transient problems, 
it Is In the latter case in which the computer has been used more frequently. 

The analog computer has been found useful in the solution of ordinary differential equations 
with constant coefficients such as the equation 


a*, - 3 -t + a —+ . .. + a —— 
1 dt 2 dfc 2 n dt n 


( 1 . 1 . 1 ) 


where the coefficients a^,...,a are constant, and where F is a forcing function which may be 
either zero or some arbitrary function of t. In fact, it Is not much harder to generalize the 


forcing function F to the form 


- f(t & 


( 1 . 1 . 2 ) 


That is, the right hand side of the equation may be a function of the independent variable t 
and one or more of the derivatives. 

The computer may also be used for solving linear differential equations with variable 


coefficients such as 


f i (t) $+... + f n( t) A= P(t), 

dt 


(1.1.3) 


where the functions f^(t) are functions of time rather than constants. 

One of the most useful applications of the analog computer is in the solution of non-linear 
equations, which in general do not lend themselves to analytic solution. As we recall, a non¬ 
linear differential equation Is one in which the dependent variable or one or more of Its der¬ 
ivatives appear raised to a power other than unity or as the argument of a function. Examples 
of non-linear ordinary differential equations are: 
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e(i-* 2 )§f+ X = o 
m + c S + K i 3ln y = k 2 

& + x c = 0. 

dtr 


Certain types of partial differential equations, linear and non-linear, can also be 
solved with the electronic analog computer. The most important are: 


a) 


b) 


Parabolic equations, such as the ’’diffusion equation” 


v 2 <j> 


d0 


fr + k 2^ + * 3 * 


which is important in the study of heat transfer. 
Hyperbolic equations, such as the wave equation 


V 0 


= k 




1 at 2 


+ k„ 


7)0 


+ k. 


2 & t ' 

which governs many physical systems involving the propagation of waves. 
Elliptic equations such as Poisson’s equation 


a 2 0 


a 2 0 


c> X 2 a y 2 


= k 


can be solved more conveniently with less expensive passive analog networks and are seldom 
treated with electronic analog computers. 


1.2 The Components of an Analog Computer 
Suppose that the equation 

a 0 y + If " a 2 = 0 

is to be solved with an electronic analog computer. By observing the operations required in 
this equation, we conclude that components are needed to: 

a) multiply a variable or a function times a constant, as in 

d!y , 

a n y , a 0 —% , etc., 

0 d dt^ 

b) generate the derivatives of a variable, 

^(t), 

of variables and functions 
the equation. 


c) 

d) 

e) 


generate arbitrary functions such as 
multiply two functions such as f^(t) 
and 

perform the addition and subtraction 
as required by the left-hand side of 


- 2 - 
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In an electronic analog computer these basic operations are performed on voltages by- 
elec tronic ’’black boxes” which accept voltages at their input terminals, operate on these 


voltages continuously, and produce voltages at their output terminals that are some function 
of the input voltages. 

The analogy involved in the use of this computer is the analogy which the user must set 
up between the physical variables such as displacements, velocities, angles, pressures and the 
various voltages present in the computer. The procedure of setting up this analogy is very 
simple, and will be explained in Chapter 2. Before this is done, however, we must study the 
individual computer components which perform the required fundamental operations. 

1 .3 Multiplication by a Constant 

One of the simplest operations which may be performed on a voltage is multiplying it by 
a constant. If the constant is less than 1, the operation may be performed with a potentiometer 
or attenuator (sometimes called a "pot”) as 3hown in Pig. 1.3,1 

The circuit diagram in Pig. 1.3.l£) shows that if a voltage E 1 (t) is applied at the input 
terminal, then any fraction K of this voltage may be transmitted to the output terminal by 
moving the slider up or down. The fraction K is always less than 1, and may be indicated by a 
calibrated dial which is connected to the slider. A diagram of the dial is shown in Fig. 1.3.3(b). 

It is more convenient to think of the potentiometer in terms of its block diagram shown in 
Pig. 1.3.3(c). This is a one-line diagram (we do not show the third terminal which is connected 
to the machine reference ground) which conveys the idea of the operation performed without 
any reference to circuit details. 
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The potentiometers used in most analog computers have a resistance element which is a hel±x 
of ten turns. The turn counter in the dial displays the turn number in a small window. This 
turn number corresponds to the tenth digits of K. The face of the dial is divided into ten 
parts which represent the hundredth digits of K, and each of these is divided into ten parts 
representing the thousandth digits of K. In this manner, the constant K may be set with an 
accuracy of one part in a thousand. Pig. 1.5.3(d) shows the block diagram of the "pot" corres¬ 
ponding to the dial setting shown in Fig. 1 . 5 . lib;. 

The input and output voltages E 1 (t) and E 2 (t) are in general functions of time. In the 
condensed diagrams of Pigs, 1.3.1c and d the symbol indicating time dependence (t) is dropped 
for simplicity, but it is always implied. The same thing will be done in some of the block 
diagrams in these notes. 

The potentiometer is a passive device and can only provide attenuation, that is, the 
constant K can never be greater than 1. An active device which provides both attenuation and 
amplification is the ^onst^an^mult^ipli^er block of Fig. 1.3.2. Figure 1.3.2a shows how this 
block is implemented from an .amplifier, usually 

•fP 1 



* - ^ E,(t) - - K £,(t) 


a) 



Fig. 1.3.2 



h 


called an "operational amplifier," and two resistances R Q and R^. The operational equation 
to the right of the diagram indicates that the output voltage is equal to the negative of the 
input voltage times a constant equal to the ratio Rq/R^. This ratio may be called the gain 
of the block, and it may be adjusted to values greater than as well as less than 1. If this 
gain is K, the one-line diagrams of Fig. 1.3*2b and c are more convenient to use and are less 
circuit oriented. The notation of Fig. 1.3.2c will be used in these notes. 
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It should be noticed that the circle and the triangle of Fig. 1.3.2c form a symbol and 
should not be confused with the symbol of an operational amplifier in series with a potentiometer. 
The figure below shows the difference between these symbols. 



Fig. 1.3.3 shows four examples of the use of the "constant multiplier." In Fig. 1.3*3a 
a constant voltage of 5 machine units (5 volts) is multiplied by a gain of 3, and reversed in 
sign, to produce a constant output of -15 machine units. In Fig. 1.3.3b the input voltage is 
negative and the gain is less than unity. In Fig. 1.3.3c the time-varying function F(t) is 
multiplied by -a, and in Fig. 1.3.3d the analog block is used to change the sign of the 



Fig. 1.3.3 

voltage E. This last block occurs so often, that it has been given the special name of slgn^ 

£hanger. 

The reader who does not have a strong background in electronic circuits should not be 
confused with words such as "amplifier," "gain," "resistance," etc. The use of a modern analog 
computer does not require such a background since we can set up our problems in the form of 
logical diagrams which are circuit independent. After some practice, anyone can learn to inter¬ 
connect the resistances, capacitances, and amplifiers to implement the diagrams. The most 
important phase of the problem-solving process is the construction of the logic diagram. The 
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actual wiring of the computer components is the job of a technician, and many installat±ons 
may provide this service to the user. 

One thing must be kept in mind, however, even when we are programming at the logic dia¬ 
gram level, and this is that all analog components have certain restrictions on their operating 
ranges which must be satisfied. For example, there are limits on the size of the resistances 
R q and which may be used in connection with a particular amplifier design. There are also 
limits on the minimum and maximum gains of a block and on the minimum and maximum value of the 
input and output voltages. 

Common ranges of resistance, gain and voltage are .01 £ R 10 (megohms), 1/50 * Gain ^ 50, 

and -100+100 volts. The reader should consult the computer manual for the actual ranges 

in the machine at his disposal. 

Notice that the common values of resistance are given in units of megohms (a million ohms 
or 10^ ohms). This makes the megohm a convenient unit for analog circuits. The reader will 
often see the ohmic value given in "megs," so that a sign-changer may show the values of 
R q = R-j_ = 1 meg. 

As seen in Fig. 1.5.2, the multiplying constant K has to be adjusted by selecting the 
proper ratio of two resistances Rq and R^. Whenever two such resistances are not available, 
a potentiometer may be used in order to adjust the gain to the proper value. 

Example 1 . 5.1 

A value of K = 2.56 is needed, but two resistances which provide this ratio are not 
available. Other gains such as .01, .1, 1, 2, 5, 10, etc., are available. The procedure is to 

divide K by. one of the available gains greater than K such that a gain K f less than 1 will 

result. The gain K' may then be set in a potentiometer. Two ways of doing this are shown 
in Fig. 1.5.5 



Fig. 1.5.5 
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1.4 Analog Addition 

The addition of two machine voltages E^(t) and Eg(t) is accomplished by the adder block 
or simply the adder or summer^ The diagram of a two-input adder is given in Fig. 1.4.1a. The 
one-line diagram which will be used in these notes is given in 1.4.1b, and the block diagram of 
an n-input adder is shown in 1.4.1c. 


«> 

K, 

:-y'V'- 


e.u) 

R*. 


4 




Ro 

\AA 


AMP 



e„u>-- -(^ ejo) 


a) 



pig. 1 . 4.1 

Notice that the adder multiplies each input variable by its corresponding gain, adds these 
products and changes the sign of the sum. Fig. 1.4.2 shows two examples of the addition of 
analog variables. 



Fig. 1.4.2 
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Subtraction of analog variables may be done by changing the sign of the subtrahend with 
a sign-changer as shown in Fig. 1.4.3. 



- ( G, A - G x 8) 


1.5 Integration 

Before we discuss an analog block which may be used for integration, it is worthwhile 
to review some basic concepts of integration. First, we remember that the symbol 

r b 

\ f(x)dx (1.5.1) 

' a 

is called the defini_te_integral^ of f(x). The function f(x) is called the integrand and the 
numbers a,b are called limits_o£ Integration. Now, if f(x) is continuous in the interval 


a£x<b, and if f(x) is the derivative of another function called F(x) , 

rb . b 


f(x)dx = F(x) 


= F(b)-F(a). 


( 1 . 5 . 2 ) 


Corresponding functions f(x) and F(x) may be found in tables of inde£inite_integrals_ 
of the form 

F(x) = ^f(x)dx 

Notice that the indefinite integral symbol does not have any limits of integration. 

For example, if f(x) = cos(t) we recall from memory, or find in a table, that cos(t) is the 
derivative of sin(t). we may then write 


r b b 

\ cos(t)dt = sin(t) = sin(b)-sin(a) 


( 1 . 5 . 5 ) 


In analog computer work the most common integrals are time dependent, such as 

t it 


( f (t)dt = F(t) 
' o 


F(t) - F(o). 


(1.5.4) 


The F(o) term may be called the initial ^ndjttion of the function represented by the 
indefinite integral. 

Now we shall describe the analog block which may be used for the evaluation of integrals. 

Integration of a machine voltage with respect to time is done with the integration block 
or integrator shown in Fig. 1.5.1a. The integrator consists of an operational amplifier, a 
resistance, and a capacitance. If the voltage across the capacitor is zero (the capacitor is 
discharged) when the input voltage E-^ is applied, the output voltage will be equal to a constant 
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times the definite integral of the input voltage with respect to time. The constant is equal 
to the -1/RC ratio where R is in ohms and C is in farads, and is called the gain of the Integrator. 
Fig. 1.5.1b shows the block diagram of the l.nte£ratoi\_ In order to specify that the capacitor 
Is discharged when the integration process begins, an oval with a 0 In It is shown below the 
triangle. 




-k( e. 

10 


dt 


a ) Fig. 1.5.1 

For example: suppose that in Fig. 1.5*1 E. (t) = Sin(t) and k = 1. We obtain 

^ T . 4- 


*‘> - ■ c 


sin(t)dt = - 


] 

o J 


(—cos(t)) 

[-cos(t) + cos(O)] = oos(t) -1 
If the capacitor is not discharged at the beginning of the integration process, the voltage 
across the capacitor will be superimposed as a constant on the output of the integrator. 

For example, if the voltage across the capacitor Is 10 volts as shown in Fig. 1.5*2, the output 


sin it) 


J U/--L T XU • lw 

_-11+_ 

1 

1 IT 


*■— vW- 

■) 




- COS it) - l +10 


Fig. 1.5.2 


On the other hand, if the voltage across the capacitor Is reversed, as shown in Fig. 1.5*3> the 
output will be equal to cos(t)-l - 10. 


tW 


SINlt) 



E a 5 cosit) -1 -10 


Fig. 1.5.3 
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The action of the capacitor voltage is shown symbolically in Fig. 1.5.4 



E/t)eJt 


-V 


In this figure the capacitor voltage is indicated inside the oval. This voltage is applied to 
the capacitor by means of a switch. At the beginning of the integration process (t = o) , the 
switch is closed and the capacitor is charged to the voltage in the oval. This appears as a 
constant term added to the definite integral. As long as the switch is closed the integrator 
output voltage is held constant and equal to the voltage across the capacitor. For this reason 
the voltage in the oval is called the initial condition voltage. In some computers the Initial 
condition voltage is applied to an amplifier input called the initial condition terminal. This 
initial condition input is similar to the other amplifier inputs except that a switch in series 
with it may be used to connect or disconnect this input from the amplifier. In this case the 
amplifier reverses the sign of the Initial condition voltage so that a negative voltage -V 
will appear as a positive constant +V at the output^ as shown in Fig. 1.5*5 




Fig. 1.5.5 


SM)Jt ■<- V 


Rather than using the convention of Fig. 1.5.5, which is computer oriented, in these 
notes we will write in the oval the actual initial condition value as in Fig. 1.5.4. Users 
of computers which produce a sign reversal should take this fact into account or may prefer to 
use the diagram of Fig. 1.5*5. 

A fact that will help us later when we start using the integrator is that if the initial 
condition voltage is equal to the factor-F(o) of equation (1.5.4), then the integrator output 
is simply the indefinite Integral of the input voltage . This is shown in Fig. 1 . 5 .6 
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When using the integrator it must be remembered that the indefinite integral will appear 
as the output only if the initial condition voltage is adjusted to the value of F(o). Other¬ 
wise, the expression representing the output voltage must be found by evaluating the definite 
integral. Example 1.5.2 showed two cases in which the indefinite integral resulted because the 
initial condition voltages were set properly. If, on the other hand, the initial condition 

voltage in Example 1.5.2a is set equal to-10 volts, then the output is given by 

f t 

E Q (t) = - J sin t dt -10 = cost - 1 - 10 

> '° 

= cost - 11 

The examples of Fig. 1.5.7 show what the output of an integrator will be when one 
machine unit (1 volt) of different polarities is applied to the input. Notice that only in the 
top example may we use the indefinite integral for the evaluation of the output. 



-ut +0 



E. 


+ <T 

t+jT 


E. - -l f -|dt -S 

'O 

: t -S’ 

Fig. 1.5.7 
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Pig. 1.5.8 shows an example involving a time varying input voltage. 




= -50 cos t 
= 50 sin t 


(real time in seconds and angle in 
degrees are analogous in this example) 


Fig. 1.5.8 


In most computers the l/RC ratio is adjusted by varying the resistance rather than the 
capacitance. This is done because it is easier to measure resistances accurately than capa¬ 
citances. The most common capacitor used has a capacitance of 1 microfarad (1x10”^ farads) 
so that a ratio of unity may be obtained by using a 1 megohm (1x10^ ohms) resistor as seen in 
Pig. 1.5.9. 



^ (l oil** )(I * lO c f-arads ) 



Fig. 1.5.9 


Other ratios may be obtained by increasing or decreasing the value of R or by using 
potentiometers as was done in Example 1.5.1. 

Example 1.5.1 

Give the block diagram of an integrator with a gain of 1 . 58 . Solution: If we assume 
that integer gains of 1, 2, and 5 are available, the diagrams of Pig. 1.5.10 are equivalent 
solutions. 
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Pig. 1.5.10 


1.6 The Summer-Integrator 


The weighted sum of the integrals with respect to time of several time-varying input 
voltages may be obtained by the summer-integrator or adder-integrator in Pig. 1.6.1. 







Generated on 2015-10-13 23:48 GMT / http://hdl.handle.net/2027/mdp.39015040284898 
Public Domain, Google-digitized / http://www.hathitrust. 0 rg/access_use#pd-g 00 gle 


Examples of the use of this block are given in Pig. 1.6.2. 



1. 7 Function Generators and Multipliers 

In paragraph 1.2 we discussed the need for devices which generate arbitrary functions of 
one or more variables and devices which produce the instantaneous product of two or more variables 
which are changing as functions of time. Devices called function generators and function multi¬ 
pliers are used for these two operations. 

Function Generators : Fig. 1.7.1a shows the general 

block diagram of a function generator. In its most general form this device accepts n time- 
varying voltages v^, v 2 , * ,,,v n and P rociuces an output voltage which is a prescribed function of 
the input voltages. 

In some cases a function generator may consist of a few electronic diodes and resistances. 

In general, however, a function generator Is a complicated piece of equipment containing many 
components. In commercial models the function may be set or changed by various means, such as 
by setting dials, inserting templets with the same shape as the graph of the function, etc. 

The most common types of function generators are those which produce a function of a 
single variable as shown in Fig. 1.7.lb. and c. Fig. 1.7.1c results when time is the argument 
of the function, so that the generator needs no input voltage, since time may be supplied 
internally by integrating a constant voltage as was done in Fig. 1.5.7. 



Fig. 1.7.1 
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Common examples of function generators of the type shown in Pig. 1.7.1c are the sine- 
wave, square-wave, and triangular wave oscillators available in any electronics laboratory. 

Fig. 1.7.2 shows two examples of the generation of functions. The reader should study 
these examples carefully. 



Fig. 1.7.2 

Function Multipliers : The general block diagram of a function multiplier is shown in Fig. 1.7.3a. 
Each input variable is in general a continuous voltage changing with time. The multiplier 
develops an output voltage equal to a constant times the instantaneous product of the input 
voltages. The multiplying constant K may be either positive or negative. In these notes we 
will show the constant inside the block. 

Fig. 1.7.3b shows a common type of function multiplier. This is a special case of the 
general multiplier which allows the multiplication of two variables only, and has a scale factor 
of -.01. 



Fig. 1.7.3 


Another common type of multiplier might be represented by the block diagram of figure 1.7.3c, 
in which a variable Z may be multiplied by each of n other variables y^,...,y n . This multiplier, 
usually of the "servo" type, then provides n outputs which are proportional to the product 
of Z with each of the n inputs. 
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An example of the use of the two-input multiplier is given in Pig. 1.7.4. Two input 
voltages v^ = 10sin9 and v^ = 10© are multiplied instantaneously to produce an output equal 
to -.01 (10 sin 9)(10 0) = -0 sin 9. 


V,- IOSfA/(d)-* \ 

V x = \OG -* - o, 


© S IN 10) 


Fig. 1.7.4 



Fig. 1.7.5 

Notice that in the function v^ = 10sin9 the analogy between radians and volts is implied. When 
the angle 9 is equal to tt/ 2 radians its analog is equal to 1.5708 volts and the sine of the 
angle is equal to 1, so that the voltage at the v^ terminal would be 10 volts. If an ordinary 
d-c voltmeter is connected between ground and this terminal, it would read 10 volts at the time 
when 9 is 1.5708 volts. 

Many methods have been developed for obtaining function multiplication. Perhaps the 
most popular methods are provided by the servo-multiplier and the square-law electronic 
multiplier. The most common values for the constant K in Fig. 1.7.1 are + 1 and + .01, depend¬ 
ing on the particular multiplier. The accuracy with which multiplication is performed depends 
on the frequency of the input voltages. Static accuracies of .1% are common and accuracies 
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of . 01$ or better may be obtained in more expensive equipment. The dynamic accuracy deterior¬ 
ates as the frequency increases, and it is better in the electronic multiplier than in the 
servo-multiplier. A detailed explanation of the different methods for obtaining function 
multiplication is not within the scope of these notes. 

The reader might notice that the function generator is the most general analog block, 
since all of the blocks discussed are special cases of it. For example, an integrator is a 
function generator which accepts a voltage and generates its time integral, an adder accepts 
several input voltages and develops their sum, a multiplier is similar to an adder except that 
it develops the product of the input variables, etc. As a matter of fact, we showed in Fig. 
1.5.7 how the functions t, t+5, t-5 may be generated with an integrator. 

A Multiplier used for Division 

A function multiplier may be used for division by connecting it to an amplifier as shown 
in Fig. 1.7.6a. Notice that the amplifier serves as an adder whose output voltage is fed back 
to the multiplier. Fig. 1.7.6b gives the block diagram of a divider. 





EM) 

EM 


E, 

b) 

pig. 1 . 7.6 

1.8 Other Analog Blocks 

There are many analog computer components which we have not mentioned here. The blocks 
which have been presented, however, are enough for solving a large percentage of the problems 
which are tackled by analog techniques. 

The output devices, that is, meters, recorders, etc. which are needed to obtain the final 
answers are explained in Chapter 3, where a sample computer is described. 
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Exercise 1.1 


DRILL EXERCISE ON ANALOG BLOCKS: 


a. 



b. 



-asv 
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1 . 




output voltage (as a 
function of time) = 


output voltage = 


output voltage = 


output voltage = 
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missing gains and 
Initial condition; 


M}/ tolling gains 




( & + ) ■+ i from the functions ® s f> 


1,1 Generate tfte function 
a single analog block. 


by using 


1,5 Notice that 


Use this idea- tci'\generate/^i^ctiaa multiplication from two Identical function generators 
••which. generate:-th& ; squaring..'function plus any .additional analog blocks needed to form 
the terms to te squared v 

1.} Draw the block dlagrson of an analog circuit for generating the function -t < in{B)> 

assuming that the function v 0 =- in() may be Implemented With a function generator* 
which accepts the voltage and produces the voltage v Q . 

l r 4 Dray* the block diagram of an analog circuit which will accept 5 variables x 7> ,,, 

4 * , /£* f ' • 5 

and generate jH v ~ V i-«» / . Generalize the block diagram to accept n variables, 

1.5 A certain analog computer has operational amplifiers which are pre-wired as shown 

In the drawing* The amplifier has a selector switch 5 Which converts It into ah adder 
if the switch is In position A, a summing integrator if in positions XI or 1,1, The 
foilowing jacks are available for connections i the seven Input lacks, four common 
output jacks, a Junction jack, and an initial condition Jack, Inputs which are not needed- 
may be left disconnected. Thus, this amplifier may fee used as a urdveraal analog blocki. 
Constant 1 ’megohm resistors'are available; and. may be coruiected to the junction, jack 
when more Inputs are needed. 
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1.5 (cont'd) 



Show how the following blocks may be implemented: 

a. sign changer 

b. constant multiplier (How many different constants can you 
get without the use of attenuators?). 

c. three input adder with gains 1,2,5 (Hint: resistances may 
be paralleled or attenuators may be used). 

d. an integrator with unity gain. 

e. a summing integrator with gains 1,.5*2 
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Chapter 2 ANALOG SOLUTION OP EQUATIONS 


2.1 Introduction 

In Chapter 1 we discussed the basic operations which may be performed on voltages by the 
components of analog computer Now we shall show how an analogy may be set up between the voltage 
in the computer and the variables in certain types of problems. 

In the discussion which follows we shall prefer to use the block diagram of the analog 
components rather than the circuit diagram. By doing this, the diagrams will be useful regard¬ 
less of the type of analog computer on which the problem is solved. 

2.2 Linear Equations 

The simplest equation that may be solved in an analog computer is a linear equation. For 
example, the equation 

y = mx + b, (2.2.1) 

which is recognized as the equation of a straight line with intercept b and slope m, may be 
generated by means of the analog block in Pig. 2.2.1. 



Fig. 2.2.1 

The equation may be solved by assuming that two voltages X and B are analogous to the 
mathematical variables x and b. If these voltages are fed to the two Inputs of an adder as 
shown in Pig. 2.2.1 and the gains are adjusted as shown in the figure, the output of the adder 
will be equal to the negative of the sum of the inputs times their repective gains, which ac¬ 
cording to the given linear equation, is equal to -Y. The output of the sign changer Is then 
a voltage analogous to the mathematical variable y. The voltages X, B, and Y are the machine 
variables, and they are analogous to the problem variable^ x, b, and y. 

Other computer connections which yield the solution of (2.2.1) are given In Pig. 2.2.2. 

The reader should study the diagrams and prove to himself that they only differ from Fig. 2.2.1 
in the way in which the machine variables are defined and distributed. For example, in 2.2.2a 
the problem variable b is made analogous to the gain B rather than to a voltage, and in 
2.2.2d negative voltages are used for X and B in order to eliminate the sign changer. This is 
easily done, since analog computers have a voltage supply which provides both polarities of 
voltage. 
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Preference for one of the diagrams in Fig. 2.2.2 depends on whether it is easier to change 
tine magnitude of a voltage or the magnitude of a gain. This depends on the computer which is 
being used. 




Fig. 2.2.2 


The problem variables x, y, m, and b are given in units which depend on the physical 
problem for which the equation is a mathematical model. For example, one form of the linear 
equation which we remember from elementary physics is , 

y = Vt + y 0 (2.2.2) 
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which gives the total distance y traveled during the time t by an object with a constant velocity 
v and an initial distance y Q from a reference point. If y and y Q are in miles, v in miles per 
hour and t in hours, the analogy between the physical and problem variables in Pig. 2.2.3 will be 

voltage T is analogous to hrs, of time t 
voltage Y q is analogous to miles y Q 
voltage Y is analogous to miles y 

non-dimensional gain v is analogous to miles/hr v 



Now we may ask: What voltages T, Y, and Y q are equivalent to hrs. of t, and miles of y and y Q ? 
Obviously, we must set up proportionality constants between each machine variable and each 
problem variable so that the machine voltages will always be within the range for which the 
analog blocks were designed. If for a particular problem t is equal to 2000 hrs, it is 
obvious that a one-to-one correspondence between t and T might damage an analog block which Is 
designed to operate at a maximum voltage of + 100 volts, since this would necessitate feeding 
2000 volts to one of the inputs of amplifier 1 in Fig. 2.2.3. 

Let us investigate in more detail how the relationships between problem variables and 
machine variables may be found. We will do this with a simple example. 

Example 2.2.1 

Given: v = 2miles/ hour 

y 0 = 100 miles 

Find: y for 0 t * 200 hrs. 

Solution: Equation (2.2.2) becomes 

y = 2t + 100 (2.2.3) 

The maximum value of y is 


'i = 2 t + 100 
max max 

= 2 x 200 + 100 = 500 miles 


(2.2.4) 


If this equation is to be solved by analog techniques, the machine variables Y, T, and Y q must 
never exceed the maximum operation voltages or gains of the analog blocks. If a one-to-one 
correspondence is set up between machine and problem variables, and if the diagram of Fig. 2.2.3 
is used, the output of amplifier 2 would be required to deliver 500 volts and one of the inputs 
of amplifier 1 (T) would require a maximum voltage of 200 volts. Assuming that our computer 
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does not handle voltages above + 100 volts, we must use a set of appropriate scal£ factors, rather 
than a one-to-one relationship between problem and machine variables. For example, since y 
reaches a maximum value of 500, each machine unit of Y (volts) may be made equivalent to 5 
problem units of y (miles), each machine unit of T (volts) equivalent to 2 problem units of 
t (hrs.) This is written 
( 2T = t, 

. 

5Y = y, (2.2.5) 

and corresponds to the process of changing variables in mathematics. 

When these relationships are substituted into (2.2.5) we obtain 

5Y = 2'(2T) + 100 

and 

0^2T<200, (2.2.6) 

which after simplification, yield the machine equations 

( 4 

Y = -jr • T + 20 

o £ T * 100, (2.2.7) 

and the circuit given in Fig. 2.2.4. 



Fig. 2.2.4 

We must keep in mind that in order to interpret the analog results, the machine units 
of volts must be changed into problem units by applying the scale factors in reverse. For 
example, when T = 50 volts, the output of amplifier 2 in Fig. 2.2.4 will be given by (2.2.7) 
and it will be equal to 60 volts. The machine variables are given by 
t =2T=2x50=100hrs. 
y=5Y=5 x 60=500 miles 

Another way to interpret the change in scale is to draw the block diagram shown"in 
Fig. 2.2.6. Here the relationships of (2.2.5) are used to label the block diagram in terms 
of the original problem variables. 



Fig. 2.2.6 
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Notice that in this diagram the actual problem variables t, y, and y Q are shown rather than the 
scaled machine variables T, Y, and Y Q . When this approach is followed, the analog circuit is 
adjusted in such a way that the outputs of the blocks will yield fractions of the problem 
variables rather than the variables themselves. In this way we avoid the overloading of the 
individual blocks. 

Notice, for example, that amplifier No. 1 would yield 

'V = ' ( \ 'I + *§) (2.2.8) 

or 

y = 2 t + y 0 ( 2 . 2 . 9 ) 

which is indeed the given linear equation. When y = y _ = 500 miles, the output of amplifier 

max 

1 will be y/5 = 500/5 = 100 volts. If this output is recorded, the recorded result should be 
labeled y /5 and interpreted accordingly. 

The diagram of Fig. 2.2.6 may be obtained by multiplying and dividing each variable to be 
scaled by the same scale factor so that the equation is not affected. In the case of equation 
( 2 . 2 . 5 ) we would like to use a scale factor of 5 for y and a factor of 2 for t, so that the 
equation may be written 


Mi) - 2 - Ml) 

The term is then Isolated to give 

( 9-1 (!)♦*>■ 


+ 100 


( 2 . 2 . 10 ) 

( 2 . 2 . 11 ) 


which is the equation represented by Fig. 2.2.6. 

Another example will further illustrate this procedure. 

Example 2.2.2 

Draw the block diagram of an analog circuit for finding the value of P as a function 
of Q and R given by the equation 

P = 400Q + .02R - 500, 

for 

04Q^2 
0$ Ft i 10^ 

1) The maximum values of Q and R are specified by (2.2.15) and (2.2.14), 

minimum values of P are pe found to be 

P . = 400-0 + .02-0 - 500 = -500 

min 

P = 400-2 + .02 x 10 4 - 500 = +500 
max 

2) The three variables must be scaled. This may be done in many different ways. One pos¬ 
sibility is to draw an analog circuit which will display the variables 50Q, r/100 , and P/5. If 

this is done, the maximum voltages represented by these variables will never exceed |+100 volt s 
The equation is rewritten in the form 


( 2 . 2 . 12 ) 

(2.2.15) 

(2.2.14) 

The maximum and 
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400* 


1 

50 


(50Q) + .02 •10' 




500 , 


which yields after simplification 

(!)-§• (5CX » + f(i?) - 100 

3) The diagram of Fig. 2.2.7 is then drawn. 


(2.2.15) 



4 

Notice that the maximum positive value of P is obtained when Q = 2 and R = 10 . In terms 

A 

of our scaled diagram, this occurs when (50Q) = 100 volts and = 100 volts. The diagram of 

Fig. 2.2.8 shows this condition. 



Fig. 2.2.8 

Another Solution : 

R P 

Suppose that for some reason we wish to display 40Q, 75 ^ , and —jq . In this case 

the equations written in terms of these scaled variables become 

10 Ho) = 400 ‘ io (40Q) + - 02 - 200 (255) - 5 ° 0 > 

( TO) = ( 40+ - 4 ’(“20o) _50, 

The diagram Is shown In Pig. 2.2.9. 

wo$ —(T)_' s, n n 

^ ZqI/ 

Fig. 2.2.9 
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The diagram shown in Pig. 2.2.9 will display a maximum voltage of -50 volts when Q and 
R reach the given maximum values. 


j 2.5 Solution of Differential Equations 

The solution of differential equations with the analog computer will be explained with 
an example of a simple ordinary differential equation with constant coefficients. 

Consider the differential equation 


a 2 + a i af + V = 0 

dt 


(2.5.D 


and the Initial conditions 


y (°) = ^00 

SSL -*10 

I t- O 


( 2 . 5 . 2 ) 

(2.5.5) 


These equations may be written in "dot” notation as follows: 

a 2 y(t) + a ± y(t) + a Q y(t) =0 (2.5.4) 

y(0) = y Q0 (2.5.5) 

y(°) = y 10 . (2.5.6) 

The first step in solving the equation is to set up an analogy between the problem variables 
and the machine variables or voltages. This is done by assuming that in some way we will be 

able to obtain a computer voltage which will change with time in the same fashion as y(t) 

changes with respect to time t. In the same manner, we will assume that voltages proportional 

to the first and second time derivatives of y are available in the computer. The second step 

is to draw an analog diagram which will relate these assumed voltages in such a way that 
equation (2.5.1) and its initial conditions are satisfied. 

One of the most useful methods of finding the analog circuits of differential equations 
is the "highest derivative method". The method consists of solving for the highest derivative, 
in this case y (t ), to obtain 


y(t) 



y(t) + 



y(t) 


( 2 . 5 . 7 ) 


and to assume that a machine voltage analogous to this derivative appears as the output of 
an adder which satisfies equation (2.5.7) as shown in Pig. 2.5.1 



Fig. 2.5.1 
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We then realize that voltages proportional to y(t) and y(t) are needed in order to 
satisfy the inputs of the adder. These voltages are obtained by repeated integration of the 
highest derivative as shown in Fig. 2.3.2 



Fig. 2.3.2 

Due to the change of sign produced by the analog blocks, the first derivative y(t) appears 
with a negative sign. 

A sign-changer may be used as shown in Fig. 2.3.3 to generate the input needed by amplifier 1. 
This figure also shows that since the outputs of amplifiers 3 and 4 represent the desired inputs 
to amplifier 1, actual physical connections are provided in order to complete the circuit. 

Fig. 2.3.3 also shows how the Initial conditions of y(t) and y(t) are applied in order to complete 
the solution of the problem. 



Notice that the initial condition of amplifier 2 is not the initial value of y(t), but rather 
the initial value of -y(t), which is equal to -y 1Q . 

Operating Procedure 

After the logical diagram of Fig. 2.3.3 is drawn, the solution of the differential equation 
on the analog computer is obtained as follows: 

1) The gains of the analog blocks are adjusted by choosing resistors and capacitors of the 
proper values or by using attenuators. 

2) The blocks are interconnected by means of pluggable wires. 

3) The initial condition voltages-y 1Q and y QQ are adjusted and applied to the initial condition 
terminals of the Integrators. If the signs of these voltages are reversed by the amplifiers, this 
must be taken into account. 
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4) The initial condition switches are closed momentarily to charge the integrating capacitors. 
When this is done, the outputs of the blocks are 

Output of amplifier 1 = - ^ y 10 + ^ y 00 ) 

Output of amplifier 2 = - y 1Q 
Output of amplifier 3 = Yqq 

5) The initial condition switches are opened, permitting the integrators to start operating 
as such. The instant of time when the switches are opened is analogous to t =0. The voltages 
y(t), -y(t), and y(t) constitute the solution to the problem, and may be recorded on a voltage 
recording device such as a strip-chart recorder. 

The analog computer yields not only the dependent variable y(t), but also its derivatives. 
The solution may be repeated as many times as desired by returning to Step 4. When this is done, 
it is possible to make changes in the parameters and the initial conditions by merely changing 
the gain of the blocks or the values of the initial condition voltages. A complete set of 
solutions may be run in a few minutes, since changing gains or voltage magnitudes is as easy 
as turning an attenuator dial. 

2.4 The Non-Uniqueness of the Analog Circuits 

If instead of solving for the highest derivative of (2.3.4) we would have chosen to solve 
for a 2 y(t), the equation 

a 2 y(t) = - a x y(t) - a Q y(t) 

would have resulted. It may be easily verified by following the same line of reasoning used 
in the previous section that this equation is solved by the diagram shown in Fig. 2.4.1. 



Other diagrams that may be used are given in Fig. 2.4.2. The reader should study these diagrams 
carefully and try to find other variations. 


30 


Digitized by 


Gck igle 


Original from 

UNIVERSITY OF MICHIGAN 



Generated on 2015-10-13 23:50 GMT / http://hdl.handle.net/2027/mdp.39015040284898 
Public Domain, Google-digitized / http://www.hathitrust. 0 rg/access_use#pd-g 00 gle 




b) 

Pig. 2.4.2 


The diagrams of Figs. 2.4.1 and 2.4.2 point out the fact that analog computer circuits are not 
unique, and that any one equation may be solved by many different circuits. A particular solution, 
however, may be more convenient in a particular case. For example, in Fig. 2.4.2a the gains of 
amplifier 1 are independent of the coefficient a^. In some particular case this circuit may 
be more desirable than any other of the many possible circuits. 

2.5 Scaling a Differential Equation 

In section 2.2 we saw how a linear equation was scaled by choosing a proper set of display 
variables which guaranteed that the outputs of the amplifiers would operate within their useful 
range. 

Scaling a differential equation is similar to scaling a linear equation, except that two 
new problems have to be solved. The first problem is that in differential equations it is some¬ 
times necessary to change the scale of the independent variable, that is, to change the time 
scale, in order to operate within the frequency range of the analog devices. This problem is 
treated in paragraph 2.6. The second problem is that it is more difficult to estimate the 
maximum values of the variables in a differential equation than in a linear equation. This 
problem is treated in more detail in paragraph 2.7. 

We shall first treat the problem of magnitude scaling, that is, the problem of keeping 
the magnitude of all the analog voltages within certain limits. The idea behind magnitude 
scaling is again to find a set of display variables which will make use of the full range of 
operation of the amplifiers and other analog components. The form of the display variables 
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again depends on the maximum values of the original problem variables. We shall illustrate 
the scaling problem by giving an example. 

Example 2.5.1 

The equation 

'x(t) + 2 x(t) + 5 x(t) =0 ( 2 . 5 . 1 ) 

is to be solved with the initial conditions 
x(0) = 20 In 

x(0) = 20 In/sec. (2.5.2) 

Assume that from the knowledge of the physical system the maximum values of the variables 
have been estimated to be 

l x(t) l max = 20 in 

\x(t)l max = 200 in/sec ( 2 . 5 . 5 ) 

• . o 

\x(t)\ max = 500 in/sec 

The analogy is now set up between a set of voltages in the computer and the problem variables 
x(t), x(t), and x(t) such that voltage x(t) in the computer is equivalent to inches x(t) in 
the problem, voltage x(t) in the computer is equivalent to inches per second x(t) in the problem, 
etc. If this is the case, however, we will not be able to display x(t) and x(t), since their 
maximum values would exceed + 100 volts. At the same time the amplifier supplying x(t) will 
only reach a maximum value of 20 volts, so that its full output range will not be utilized. The 
solution is to decide on a set of display variables which will make the maximum value of each 
voltage approximately equal to + 100 volts. Such a set of display variables is: 


100 

“2S 


x(t), 


100 ;m 
SUU l ' t> ’ 


and 


100 

500 


x(t) 


5 x(t), x(t), and i x(t). 

The differential equation and its initial conditions are then written in terms of the 


display variables as follows: 

5 (^-) + 2 ■ 2 ( —) + 5 * 3 (5 x(t)) = 0 , (2.5.4) 

\ (5 X (0)) = 20, (2.5.5) 

2 (*^1) = 20. (2.5.6) 

Simplifying and isolating the highest derivative we obtain 

+ 3 + 1 (5 x( t)) = 0 , (2.5.7) 

(5 x( 0 )) = 100 , ( 2 . 5 . 8 ) 

(k2l) = io. (2.5.9) 
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The scaled highest derivative (x(t)/3) is assumed to exist at the output of amplifier 1 
or Pig. 2.5.1. This output is integrated to obtain the next lower scaled derivative (x(t)/2) 
which will appear at the output of amplifier 2. Another integration will produce the scaled 
solution (5 x( t)). 

Notice that the gains of the integrators are adjusted so that their outputs and inputs will 
correspond to the scaled variables. 



Fig. 2.5.1 

The reader should convince himself that the gains of the integrators should be 5/2 and 10 
respectively. 

It is of value now to consider how the diagram shown in Fig. 2.5.1 may be used to obtain 
a computer solution. Since the user must implement the amplifier gains by means of resistors 
and capacitors, the first step leading to a computer circuit is to draw a more complete diagram 
which may be called the circuit diagram. This is shown in Fig. 2.5.2 


I 1 I 



Fig. 2.5.2 assumes that potentiometers and resistances of .1, .5* and 1 megs, and capacit¬ 
ances of 1 microfarad are available. All the amplifiers and potentiometers should be numbered 
with the same numbers used in the computer. In this diagram it is assumed that the amplifiers 
reverse the sign of the initial condition voltages. 
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Operating Procedure 


After the circuit has been connected as shown in the figure, switches s 1 and s^ are closed 
simultaneously so that the initial conditions may be applied. The voltage at the output of each 
amplifier should then be measured to be sure that all the connections and gains have been 
properly set. In this particular example, it is easy to verify that the amplifier outputs 
should be: 


Amplifier 

1 

46.666 

volts 

11 

2 

-10.0 

volts 

m 

3 

100.0 

volts 

11 

4 

10.0 

volts 


The solution is then obtained by opening switches s 1 and s^ simultaneously, permitting 
the analog blocks to operate on voltages appearing at their inputs. Since there are no con¬ 
straints on the output voltages of these blocks^ as soon as s^ and s^ are opened they will 
change in accordance with the operational equations. The output voltage of each block is now 
measured, and if any of the four amplifier outputs exceeds the maximum values specified by the 
manufacturer, the display variables are changed and the equation is re-scaled. For example. 

If the output of amplifier 1 exceeds the maximum value specified by the manufacturer, we may 
change the scale factor of the second derivative and display x(t)/5 rather than x( t)/3. This 
would naturally require a change in the gains of amplifiers 1 and 2. 

Another difficulty that might occur is that, although none of the amplifiers is over¬ 
loaded, one or more of them may not make use of their full range of operation. For example, if 
the first derivative swings from -5 to +5 volts, we might consider re-scaling the problem so 
that a larger swing is obtained. 

If none of the amplifiers is overloaded and if all the voltages are of the proper order 
of magnitude, the solution is ready for recording. Fig. 2.5*3 shows the results of the solution 
given by the diagram of Fig. 2.5.2. 

The recorder used in this case was a stripchart recorder. Notice that each trace was 
properly labeled with the following: 

1) The name of the scaled variable. This permits the conversion between 
volts and the proper problem units. 

2) The chart speed. This permits the proper Interpretation of the time axis. 

3) The time scale. This provides information in the case of a solution which may 
have been slowed down or speeded up as explained in paragraph 2.6. 
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.The solutions obtained In an ‘analog ooci tputex' are lime dependent, that is, the dependent 
variables and' their derivatives. are .en&hglnig; ^i'tb respect to time. An important property of 
these variables -is-• theirfrequency of oscillation..,. the common u^lfca of frequency being the cycles/ 
aectrid an<X tm ra&lans/aeQdnd* 

The frequency of. oscillation la very Important because each component of an electronic 
-.analog computer ms a certain^ frequency range within which it will operate - -with; tne accuracy 
specified by the manufactureri. For example, an adder or an integrator may operate properly from 
zero frequency (consVoltages) to 500 Cycles per second, but a strip-chart recorder might 
only be abl« to operate from 0 to 30 cycles.per second. 

Consider the. computer solution of Fig. 2/6. la, We notice that the frequency of as di¬ 
lation waff bo high that it is difficult to separate one cycle of oscillation from another:. In 
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such a solution, we might want to be able to measure the values of the variable during the first 
cycle of operation. This first cycle was obtained in Pig. 2.6.1b by slowing down the solution 
so that instead of displaying the first cycle in 1 second of computer time, we display it in 
10 seconds of computer time. Since it takes longer to display the same part of a solution, we 
say the computer has been slowed-down, in this case by a factor of 10. This is similar to 
slowing down the motor in a movie projector in order to see a picture in slow motion. It 
takes longer to see the same number of frames. 




b) 


Fig. 2.6.1 

Notice that the time axis of Fig. 2.6.1b is labelled lOt to indicate that the values 
of time read from the graph are ten times larger than the actual values of problem time. The 
user must divide computer time by 10 in order to obtain problem time. If computer time is called 
T and problem time t, the relationship between them is given by 

T = lOt 

Thus, to obtain problem time t, the relation used is 

r 


t = 


10 


The solution shown in Fig. 2.6.2a, on the other hand, may be so slowly varying, that we 
might have to run the computer for five or ten minutes to get all of the values needed. This 
would be the case, for example, if we were studying some chemical process which has a slow 
rate of reaction. In this case we should speed-up the solution as shown in Fig. 2.6.2b. 



Fig. 2.6.2 
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In order to compare these figures, a reference time t^ was marked in the time axis. In 
fch.e unsealed solution shown in Fig. 2.6.2a it would take the computer 80 seconds to reach this 
reference time. In the scaled solution of Fig. 2.6.2b it would take the computer only 8 seconds 
to reach the same reference. The solution is then said to have been speeded-up by a factor 
or 10. In order to indicate this fact, the time axis of this solution is marked t/10 so that 
when the chart is read, we may realize that time must be multiplied by 10 in order to obtain 
actual problem time. 

Since many physical problems are such that the frequencies involved exceed the acceptable 
ranges in the computer, the scale of the independent variable, time, must be changed in order 
to solve these problems. 

The simplest method of changing the time scale is to proceed as follows: 

1) The equation is magnitude-scaled by using the estimated maximum values of the 
dependent variables. This was explained in 2.5. 

2) The frequencies of oscillation of the variables are estimated from knowledge 

of the physical system or they are observed after running the un-scaled solutions 
on the computer. 

5) If the frequencies of oscillation are not within the operating range of the 

computer components, the time scale is changed as follows: 

a) If the frequencies are too high , the computer solution may be slowed- 
down by a factor K if the gains of all the integrators are divided by K. 

b) If the frequencies are too low , the computer solution may be speeded-up 
by a factor K if the gains of all the integrators are multiplied by K. 

For example, if the frequencies of oscillation in a physical problem are of the order of 100 
cycles per second and the available recorder is only accurate to 10 cps, the solution may be 
slowed-down by dividing the gain of every integrator by 10. This may be done by increasing 
either the capacitance or the resistance by a factor of 10. 

If, on the other hand, the problem involves a very slow process which has an oscillating 
frequency of .01 cycles per second, the computer solution may be speeded-up by a factor of 
100 by multiplying the gains of the Integrators by 100. This may be done by decreasing either 
the capacitance or resistance of each integrator by a factor of 100. 

In some analog computers the value of the integrating capacitors may be increased or 
decreased by a factor of 10 by merely throwing a switch. This simplifies time scaling con¬ 
siderably, since the user may speed-up or slow-down a solution with the least amount of alge¬ 
braic manipulations. It should be remembered that a change in time scale by this method does not 
change the form of the variables displayed by the integrators. For example, if the output of 
an integrator is (x(t)/k), it remains the same although the input-output relationship with respect 
to time of the integrator is different. 
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2.7 Estimation of Maximum Values 


As was seen in the previous sections, scaling depends on the maximum values of the 
dependent variable and its derivatives and on the maximum value of their frequencies. 

The most successful method of scaling results when the maximum values are estimated 
from our Intimate knowledge of the sytem which we are trying to simulate in the computer. For 
example, an automotive engineer studying an automobile suspension system knows from previous 
experience the approximate maximum value of the displacement, velocity, and acceleration of an 
automobile tire. An electrical engineer familiar with an electrical network is able to make 
very good estimates of the maximum voltages and currents involved in the system. 

Whenever the scientist is confronted with a physical system with which he is not familiar 
or which has not yet been built, the problem of scaling becomes one of wits and trial and error. 

Although there is no general method for estimating the maximum values of the analog 
variables, a few methods which are successful in the case of simple linear equations have been 
worked out* Since all methods are based on some rather restricting assumptions, they may or 
may not work when applied to a particular equation. 

It is interesting to investigate the relationship between the scale factors and the 
amplifier gains. For example, the system 
a g x(t) + a^t) + a Q x(t) = 



x( 0 ) = x Q , 
x( 0 ) = Xq , 

may be scaled by choosing display variables 




/ x(t) 

l s l 


and 


'kii\ 

s 2 ) 


[2.7.1) 


where s Q , s^, s^ are 


scale factors which are usually chosen from knowledge of the maximum values of the variables. 
This yields the system of (2.7.2). 


M + !i !i 


f *(t) \ 

l S 2 ) 


+ !o !o 

3 2 a 2 


iLU 

s 0 


= 0 



f (2.7.2) 
i 

J 


The analog diagram for this equation is given in Fig. 2.7.1. If this is compared to the 
diagram of Fig. 2.3.3* which corresponds to the unsealed system of (2.7.1), we notice that each 
gain, except that of the sign changer, is a function of the scale factors s Q , s 1 , and sThe 
initial conditions are also a function of these factors. 


*For example see Analog Computer Techniques , Clarence L. Johnson, McGraw-Hill Book Co., 
New York, 1956, and 

Analog Computation, Albert S. Jackson, McGraw-Hill Book Co., New York, i 960 . 
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This suggests a very convenient method for scaling a differential equation by trial and error. 

The method may be stated as follows: 

1. Obtain the general diagram of the differential equation in terms of the scale factors. 

2. Assume values for Sq, s^, and s^ from the estimated maximum values. Use these 

values in Fig. 2.7.1 to obtain a computer circuit for the first trial. 

5. Three things may happen: 

a) No amplifier is overloaded. In this case measure the analog variables to see 
if they are of the proper magnitude for full accuracy. If they are, the trial 
and error procedure is finished. If one or more of the variables are too 
small for accuracy, change the scale factor and make another trial. 

b) One or more amplifiers are overloaded. In this case change the scale factor 
of the variables which overload an amplifier and make another trial. 

c) No amplifier is overloaded but the frequency of oscillation of the variables 
is either too high or too low. In this case we proceed as in a, except that 
the integrator gains are multiplied or divided by a constant in order to 
change the time scale. 

Example 2.7.1. 

In this example we will solve equation (2.5.1) by trial and error. Since this equation 
is a special case of (2.7.1)* the general diagram of Fig. 2.7.1 applies. 

In the first trial, we will assume that s 0 = = s 2 = When these values are written 

in Fig. 2.7.1* together with a Q =5* a i = 2, and a^ = 1, we obtain the diagram shown in 
Fig. 2.7.2. 
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If the circuit diagrammed in Fig. 2.7.2 is implemented on an analog computer, amplifier 1 
will be overloaded. The reader should verify this with the computer at his disposal. 

The fact that an amplifier is overloaded does not give us much information, since we are 
still not able to find the maximum value the amplifier is trying to reach. The trial and error 
procedure then requires that we choose a scale factor for x (t). Suppose we choose s 1 = 5. 
Substituting this value into Fig. 2.7.1 we obtain Fig. 2.7.3 



Pig. 2.7.3 

The numbers in brackets represent the maximum value of each output as measured on the computer 
circuit. Since no amplifier is overloaded, this circuit may be used for the solution. Alternately, 
these maximum values may be used to obtain coefficients which allow the computer to use more of 
the range of amplifiers 1 and 2. Such a solution may be obtained if s Q = J , s^ = ^ , and s 2 = 


2.8 Differential Equations with Forcing Functions 

If we add a forcing function to (2.2.1), we obtain 

a 2 + a i St + v = f ( t) - 

dt 


( 2 . 8 . 1 ) 


Writing this equation In "highest derivative" and dot notation, we obtain 

yU) =-jji y(t) + ^ ^+-| 2 f(t) (2.8.2) 

This is almost identical to (2.2.7), except that a factor proportional to the forcing 
function appears on the right side of the equation. The M highest derivative method" then requires 
that an additional input equal to - — f(t) be connected to the adder which produces y(t). This 
input may be generated by a function generator as shown in Fig. 2.8.1. 
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Pig. 2.8.2 shows the diagram for the equation 

a 2 + a i at + a o y = f (at) (2.8.3) 

in which the forcing function depends on the first derivative ^ . The idea may be extended 
to a forcing function which depends on one or more other variables. 



It should be remembered that the forcing function f(t) must also be speeded-up or slowed- 
down whenever the time scale in the computer is changed. For example, if the computer time 
scale is slowed-down by a factor of 10, the function generated by the function generator must 
have a frequency which is one-tenth the given problem frequency. 

2.9 Equations with Variable Coefficients 
The equation 

a 2 y(t) + a 1 y(t) + f(t)-y(t) = 0 (2.9.1) 

has a coefficient f(t) which is a function of time. If the usual procedure is followed, the 
diagram shown in Fig. 2.9.1 results. 
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Notice that since both f(t) and y(t) are time dependent, a function multiplier is required to 

obtain their product. If f(t) is an arbitrary function of time, a function generator is also 

needed to generate this function. If f(t) happens to be one of the available functions y(t), 

-y(t), or y(t), the function generator may be eliminated. If f(t) is an integral power of t, 
2 3 

such as t, t , t , etc., it may be generated by means of integrators as shown in Fig. 2.9.2. 



Fig. 2.9.2 


2.10 Simultaneous Differential Equations 

Simultaneous differential equations may be solved very conveniently by analog methods. 
The techniques used are the same as those used in the case of a single equation. For example, 
in the system 


a^U) + b u y(t) - b 12 x(t) = 0 


( 2 . 10 . 1 ) 


y(0) = y Q , 


( 2 . 10 . 2 ) 


a 22 x(t) - b gl y(t) + b 22 x(t) = 0 J 
y(o) =0 

x(0) = x Q , x(0) = 0 

we notice that the top equation in (2.10.1) is a second order equation in y, except that the 
coupling term -b 12 x(t) also appears. Similarly, the lower equation is a second order equation 
in x except for the -b 2 ^y(t) term. The method used here is to find the analog circuit for each 
equation as if the other did not exist, and then couple these individual circuits by means of 
the coupling terms. For example, if we start with the equation in y, isolating the highest 
derivative, we obtain Fig. 2.10.1. 
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Pig. 2.10.1 

If we proceed in the same way with the equation in x, the diagram shown in Fig. 2.10.2 results. 



Pig. 2.10.2 

The two diagrams are almost complete, except that -x(t) must be supplied to the y-system and 
-y(t) to the x-system. These are the coupling terms and voltages proportional to them are avail¬ 
able in their corresponding diagrams. We notice, however, that they appear with the wrong 
sign, that is, the y-system provides a voltage analogous to y(t), but the x-system needs -y(t). 
The same difficulty appears in the case of the x(t) term. This may be resolved by assuming 
that -x(t) is available at the output of amplifier 4 rather than x(t). This changes the signs 
of all the analog voltages of Fig. 2.10.2 as seen in the complete diagram of Fig. 2.10.3. 
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Pig. 2.10.3 


It should be remembered that the diagram of Pig. 2.10.3 is the unsealed diagram of the system. 
When an actual problem is solved, scaling must be considered. 
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PROBLEMS 

2.1 Draw the ’’highest derivative” logical diagram for the solution of the differential 
equation 

a Q x(t) - a 1 x(t) +‘i*(t) = 

with = ?00 


2. 2 Draw the logical diagram for 

3y(t) - 2y(t) + x(t) =0 
x(t) - 5x(t) - z(t) = 0 
2z(t) - z(t) - y(t) = 20 

with all Initial conditions equal to 0. How may we simplify the diagram if only the 
derivatives x(t), y(t), and z(t) are desired? 

2.5 Draw the circuit diagram for the logic diagram given below. Assume that resistors 

of .1, 1, and 10 megaohms and capacitors of .1, 1, and 10 microfarads are available, 
and use a maximum of two potentiometers. 




Digitized by 


Go 'gle 


Original from 

UNIVERSITY OF MICHIGAN 



Generated on 2015-10-13 23:52 GMT / http://hdl.handle.net/2027/mdp.39015040284898 
Public Domain, Google-digitized / http://www.hathitrust. 0 rg/access_use#pd-g 00 gle 


2.5 The following questions apply to the diagram of Problem 2.4 (a.) 

a. What is the output of each amplifier at t = 0 if f(0) = -500? 

b. What gains must be used in amplifiers 2 and 5 to slow down 
the solution by a factor of 5? 

c. What would happen if the function generator which produces f( t) 
is momentarily disconnected from the input of amplifier 1 at 

t = 0? 

2.6 Solve the equations 

x(t) 
y(t) 

with initial conditions 

x(0) 
x(0) 
y(o) 

2.7 Solve the equations 

x(t) + 4y( t) + 9x(t) = 0 
y(t) - 4i(t) + 9y(t) = 0 
with initial conditions 

x(0) = y(0) = 0 

y(o) = x(o) = -1 

2.8 Solve the equations 

5y(t) + 702y(t) - 4000x(t) = 0 
x(t) - 7y( t) + 800x( t) = 0 

and 

y(o) = 50 , y(o) = o 

x(0) = 2 , x(0) = 0 


+ 3x(t) - 2y(t) = 0 
- 6x( t) + 7y(t) = 0 

= y(0) = 0 
= 10 
= 30 
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Chapter 3 . The Analog Computer 


3.1 Introduction 

Just as the size of a digital computer is often measured by its memory unit, number of 
input-output devices, etc., there are two important factors in estimating the size aind com¬ 
putation power of an analog computer. These factors are: 

a) the number of operational amplifiers 

b) the number of non-linear elements such as multipliers and function generators. 

Small analog computers might have from 5 to 30 operational amplifiers and few, if any, 
non-linear elements. Large analog computers have hundreds of amplifiers and dozens of non¬ 
linear devices. The number of resistors, capacitors, potentiometers increases with the number 
of amplifiers. In some special applications the output equipment such as display and recording 
devices become quite specialized and may comprise a large percentage of the cost of the computer. 
The cost of analog computers may be as low as $1000 and as high as hundreds of thousands of 
dollars depending on the size of the installation. A good educational analog computer with 
5 to 10 amplifiers, 10 to 20 potentiometers, 2 or 3 mon-linear elements such as function gener¬ 
ators or multipliers, and a number of resistance and capacitance elements may be bought for 
about $2000-$4000. 

The schematic diagram of a small analog computer is given in Pig. 3.1.1. This computer 
does not exist, but it is very similar to most computers of its size. Most of the differences 
are in the manner in which the components are arranged. 

Because of the many variations in the equipment used in analog computers of different 
sizes and manufacturers, it is not within the scope of these notes to describe individual computer 
components in detail. It is important to the reader, however, to have a general understanding 
of what these components are. A brief description of the computer of Pig. 3.1.1 will be used 
to give the reader this understanding. After reading this description, the reader should 
consult the operating manual of the particular computer which is available to him and should 
locate and study its components parts. 
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Fig. 3.1.1 


3.2 The POTENTIOMETER PANEL 

This panel houses the ten-turn potentiometers and their corresponding dials. Each 
potentiometer is marked by a label such as PI, P2, etc., for identification. The three leads 
of each potentiometer may be brought out to plug-in jacks below it or to jacks at the central 
patch-board. Some of the pots may have one of their terminals internally grounded, that is, 
connected to the common computer ground. Then only the input and output (slider) terminals 
are brought out to plug-in jacks. In some computers it may be possible to apply the computer’s 
reference voltage, usually + 100 volts, to the input terminal of any potentiometer so that the 
attenuator constant may be set by reading the output voltage on a voltmeter as seen in Fig. 3.2.1. 

48 


Digitized by <^QuQie 


Original from 

UNIVERSITY OF MICHIGAN 




Generated on 2015-10-13 23:53 GMT / http://hdl.handle.net/2027/mdp.39015040284898 
Public Domain, Google-digitized / http://www.hathitrust. 0 rg/access_use#pd-g 00 gle 


loov 



voltage 

10c 


<f<31*1 


3.5 The AMPLIFIER PANEL 

This section houses the operational amplifiers which are needed for implementing the 
analog blocks. Each amplifier may be a self-contained unit which may be un-plugged from the 
computer for replacement or repair, or two or more amplifiers may be part of a multiple unit. 

The individual amplifiers may be of the multi-purpose type that may be used for any 
of the blocks discussed in previous paragraphs or they may be pre-wired as sign-changers, adders, 
integrators, etc. 

The front panel of a general purpose plug-in amplifier might be as shown in Pig. 5 .5.1. 

The neon light is used to warn the user that the output voltage has exceeded the maximum or 
minimum 



allowable voltage. Whenever the neon light is ON, the amplifier does not operate properly and 
it is said to be overloaded. In this case the problem must be re-scaled. 

X'The selector switch may be used to select the mode of operation of the amplifier, that is, 
its use as an adder, an integrator, etc., or it may be absent in certain designs. 

The balance adjustment knob is used in adjusting the electronic circuit so that a zero 
input voltage will produce a zero output voltage. This operation is called "balancing the 
amplifier." 

Some amplifiers are of the type shown in Pig. 5.3.2a where both polarities of the output 
are provided. This eliminates the need for sign-changers in the analog circuits at the expense 
of more equipment inside each individual block, since the amplifier in this case is equivalent to 
a regular amplifier and a sign-changer as shown in Pig. 5.3.2b . 
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Pig. 3.3.2b) 

As will be seen in Chapter 4, amplifiers of the double-output type would simplify the programming 
of analog computers if a matrix approach is used. The extra expense involved, however, has 
discouraged most manufacturers from offering computers using this concept. 

3.4 The CONTROL PANEL 

The most important function of the control panel is to provide a centralized location 
from which all the computer components are at the command of the operator. 

The four push buttoms marked BAL, INITIAL, COMP, and HOLD are connected to a number of 
Internal relays and switches which perform the following functions: When the BALANCE buttom Is 
pressed the Inputs of all the amplifiers are connected to the common ground (0 volts) by special 
networks. The operator may then adjust the amplifier output voltage to zero by turning the 
balance adjustment knob shown in Pig. 3.3.1. When the INITIAL button is pressed, the initial 
condition voltages applied to the initial condition terminals of the Integrators are sent to the 
integrating capacitors. Pushing the COMPUTE button releases the initial condition voltages 
and permits the analog blocks to operate according to their operational equations. When the 
HOLD button Is activated, all the integrating capacitors are disconnected from the integrators 
so that the capacitors will not be discharged as long as the computer is In the HOLD position. 
This is equivalent to "stopping time," since all of the voltages In the analog circuit will remain 
static if the integration process is halted by disconnecting the integrating capacitors. 

The HOLD condition permits the operator to make changes in the circuity such as changing 
potentiometer settings. This Is equivalent to changing the oroblem parameters after a solution 
has been partially run. 

In our model computer the control panel also houses a dc voltmeter and a meter scale 
switch. One side of the meter is connected to the ground reference. By means of a toggle switch 
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the other side may be connected to either an external jack or to a meter bus. A meter bus is 
nothing but a good conductor brought to the balance push button of each amplifier as shown in 
Fig. 5.4.1. 



Fig. 5.4.1 

The meter may then be used as 1) a regular meter through its external jack, 2) a balance 
meter when the computer is in its BALANCE mode, 5) a checking device when the amplifier is in 
INITIAL, or 4) a measuring device when the amplifier is in COMPUTE. In every case, except the 
first, only the output of the amplifier whose push button is being activated will be measured. 

In some computers the meter Is of the type called "digital voltmeter" which displays the 
voltages In digital form rather than with a pointer and a scale. 

5.5 The VOLTAGE SUPPLY PANEL 

The voltage suply panel houses electronically regulated power supplies which provide 
very stable voltages for use as constant voltages and initial condition voltages. In some cases 
the voltage supply panel provides only the reference voltages, usually + 100 volts, and these 
have to be scaled down with potentiometers as shown in Fig. 5*5.1. In other cases the voltage 
supply panel has built-in attenuators so that voltages of varying magnitudes and polarities 
may be obtained directly. 
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3.6 RESISTANCE AND CAPACITANCE PANEL 

This panel contains the resistors and capacitors which are used in the amplifiers. If the 
amplifiers are of the multi-purpose type which already contain pre-wired resistors and capacitors, 
the panel Is used to complement the R-C components already available. The terminals of all 
components are either brought out to plug-in jacks or to the central patch board. In the latter 
case, the patch board is clearly labeled with the value of each resistance in megohms and each 
capacitance in microfarads. 

3.7 FUNCTION GENERATOR PANEL 

As the name implies, this is where one or more function generators are located. Depending 
on the method used for the generation of functions, the panel might contain a number of knobs 
for setting the parameters of the function to be generated. 

3.8 FUNCTION MULTIPLIER PANEL 

This panel contains the function multipliers and all the necessary meters and adjustments 
for the calibration and operation of these devices. 


3.9 The PATCH BOARD 

The patch board is that part of the computer where all the computer elements are connected 
to each other by means of patch cords. The input and output terminals of all the amplifiers, 
resistors, capacitors, potentiometers, etc., are all wired internally to plug-in Jacks at the 
patchboard. The face of the board is clearly labeled to simplify the operation of the machine. 
In some computers the patch board is removable, so that it may be wired while the computer is 
being used by someone else. If several patch boards are available, several people may be 
getting ready to use the computer at the same time. It is a simple matter to Insert a pre-wired 
patch board into the computer and to run the solution in a matter of minutes. 
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Problems which are run frequently may be kept in separate patch boards ready to be 
inserted whenever they are needed. This is similar to having a library of digital computer 
programs on punched cards or magnetic tape. 

A removable plug board permits more efficient use of the computer because the machine is 
used only when the operator is ready to run his problem and it is not tied up while the components 
are being interconnected or when wiring errors are being traced. 

5.10 OUTPUT EQUIPMENT 

The most common pieces of output equipment are the strip-chart recorder, the X-Y plotter, 
and the oscilloscope. 

The strip-chart recorder is nothing but a dc voltmeter whose pointer is equipped to leave 
a trace on a moving strip of accurately calibrated paper. In some cases the pointer is similar 
to an ink pen which leaves an ink trace as it moves across the chart. In other cases the 
pointer has a heating element which burns a fine line on specially treated paper. In both cases 
the recorder will look similar to Pig. 5.10.1. 





iV\ j* 


The recorder in Pig. 5.10.1 Is a one-channel recorder which permits the recording of only one 
variable at a time. In most computers multi-channel recorders which may record several variables 
at once are available. 

The X-Y plotter is a device which has two degrees of freedom. The inking pen moves in 
the X direction (x-axis) as a function of the voltage applied between the ground terminal and 
the X input terminal. The pen moves In the Y direction as a function of the voltage applied 
between the ground terminal and the Y input terminal. In this fashion, the final location of 
the pen on an X-Y coordinate system will depend on the magnitude of the X and Y voltages. A 
drawing of an X-Y recorder is shown in Fig. 5.10.2. 
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As may be seen from the figure, the vertical bar V which carries the pen may move right and 
left by riding on the horizontal guides K 1 and H^. At the same time, the pen may move up and 

down by riding on the vertical bar. These motions are controlled by the voltages applied to 

the X and Y inputs respectively. Each input has a corresponding scale selector switch and a 
balance adjustment which positions the pen anywhere in the paper when both input voltages are 
zero. 

If a voltage is to be plotted versus real time, the arrangement of Pig. 5.10.3a may be 
used. In this case an integrator with -1 volt applied to its input generates the time function 

t which may be applied to the X-axis. The voltage to be plotted is then applied to the Y-axis. 




Pig. 3.10.3 

In other cases we may wish to plot one function of time f^(t) versus another fg(t). This 
may be done as shown in Fig. 3.10.3t . This second arrangement is not available in the strip 
chart recorder, since the x-axis is always moved at constant speed by the strip chart drive 
motor. 

An oscilloscope is a device which operates in a manner similar to the X-Y recorder, 
except that it has no mechanical parts and is all electronic. The plotting pen is now an electron 
beam which moves in the X-Y plane as a function of the voltages applied to the X and Y inputs. 
Instead of having a piece of paper for recording the answer, the oscilloscope produces a spot 
of light on a fluorescent tube similar to the tube in a TV set. As the spot moves across the 
tube, a lighted trace is produced which is visible for a fraction of a second. In this way the 
performance of the X-Y plotter is duplicated at electronic speeds. The lighted trace may be 
photographed by a special oscilloscope camera if the result is to be preserved. 
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Since the oscilloscope is all electronic, it operates at higher frequencies than the 
electromechanical X-Y plotter. For this reason, the oscilloscope is used in high speed 
analog computers of the repetitive type . These computers have special circuits which permit 
them to repeat the same solution several times each second. If the same trace is repeated on 
an oscilloscope several times a second, the observer may see the complete solution without 
having to wait until an inking pen traces across the paper. In this manner, parameters may be 
changed while the computer is displaying the solution on the oscilloscope . The effect of para¬ 
meter changes on the solution is therefore immediately apparent on this type of computer. 
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Chapter 4 MATRIX PROGRAMMING OF ANALOG COMPUTERS 


4.1 Introduction 

In previous chapters we were concerned with the use of analog blocks for the solution of 
linear and differential equations. The analog blocks were put together by using the highest 
derivative method and other methods which rely on our intuition. Although we were aware of the 
fact that several computer circuits could be found for the solution of each problem, it was 
difficult to tell if the circuit which we designed could be improved. 

In this chapter we will show a straight forward method based on the theory of matrices 
which will permit us to optimize analog circuits, to eliminate unnecessary operational ampli¬ 
fiers, and to find the relationship between equivalent circuits. In addition, the methods in 
this chapter will organize our work in such a way that our design of a good circuit will not 
depend as much on intuition. Furthermore, we will find that familiarity with the matrix 
methods in this chapter will widen our understanding of the solution of equations in general, 
and will bring to light the relationship between mathematics and the analog networks. 

Although a basic background in the theory of matrices will make some of the methods that 
follow more meaningful, the reader should not feel handicapped in reading this chapter if he 
has not been introduced to matrix algebra. The chapter will be self-contained in the sense 
that all of the necessary rules are introduced within the section in which they are needed. 
Actually, we could eliminate the word n matrix" from this chapter and consider that we are dealing 
with tables or arrays of numbers which are to be manipulated according to a number.of rules and 
procedures which yield the desired results. 

In the next few sections a few concepts of matrix algebra are reviewed very briefly. 

Only those concepts which are relevant to the treatment of matrix programming of analog computers 
are covered. A minimum number of definitions and only a few new terms are introduced. The 
reader should consult texts on matrix algebra for a formal treatment of the subject. 


4.2 The Algebra of Matrices 

In this chapter we will encounter arrays of numbers or functions which may be arranged 
into m rows anc 
by the symbol 


columns 

in the 

form of j 

a n 

a 12 • 

• • a ln 

a 21 

a 22 • 

• • a 2n 


ml 


m2 
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and may be called a matrix . A matrix of m rows and n columns Is called a m x n matrix , and the 
numbers or functions which form the matrix are called the elements of the matrix. The element 
in the ith row and jth column of the matrix is identified symbolically as a^. The elements 

a ll' 8l 22* **• 3 a mm are said to belon g to the diagonal of the matrix. 

Most of the matrices In this chapter are of two types: 

a) square matrices in which m=n, and 

b) column matrices In which n=l. 

Examples of square matrices are given in Pig. 4.2.1 



(a) (b) (c) 

Pig. 4.2.1 

The matrix shown in Pig. 4.2.1a is a 5 x 3 matrix in which every element is equal to zero. A 
matrix in which every element a^O is called a null or zero matrix, thus this matrix is a 
3x3 zero matrix. Matrix b has each diagonal element equal to unity and each non-diagonal 
element equal to zero, and is a special case of the unit matrix, in this case a 4 x 4 unit matrix. 
Matrix c is a 2 x 2 matrix in which some of the elements are expressions of the algebraic 
variables a, b, and c. 

Examples of column matrices are given in Figure 4.2.2 



Pig. 4.2.2 

The column matrix shown in Pig. 4.2.2a Is a column matrix of real constants, matrix b is a 
column matrix of variables x-^, x^, ... , x^, and matrix of c is a column matrix of functions 
of time. Column matrices are sometimes called vectors . 

The manipulation of matrices is simplified if they are given names or symbols. Matrices are 
usually named with letters of the alphabet. In these notes, boldface capital letters will be 
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reserved for naming matrices. For example, the unit matrix will be called by the name X , 
and the zero matrix by the name ©. Other matrices will be named with the letters a , B,c 
D , etc. 

The set of rules, theorems, postulates, etc., for manipulating matrices is called the 
algebra of matrices or matrix algebra. It is beyond the scope of these notes to give even a 
brief treatment of this useful and interesting subject. It will suffice to state some simple 
definitions and to review the rules of addition, subtraction, and multiplication of matrices. 
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in which 


°ij - a ij + b u 


(4.4.2) 


We may write symbolically 

C - A * ^ (4.4.5) 

According to this definition, it is obvious thatA and Bmust have the same number of rows and 
columns, i.e., both matrices must be of the same order. Similarly, the operation 

c. - A - B ( 4 . 4 . 4 ) 

is defined by the equation 


c u = a u 




(4.4.5) 


Addition of matrices which are equal leads to the concept of multiplicity, 
the addition 


For example. 



' a n 

a l2 

a l?" 


" a ll 

a i2 

a 13 


2a 

11 

2a l2 

2a 13 


a 21 

a 22 

a 25 

+ 

a 21 

a 22 

a 25 


2a 

21 

2a 22 

2a 25 


a 5i 

a 32 

a 53_ 


a 31 

a 52 

a 53. 


2a 

51 

2a 52 

2a 33_ 

suggests 

the 

notation 











" a n 

a 12 

a 13 


" a il 

a l2 

a 13 _ 



~ a u 

a 12 

a 13~ 


a 21 

a 22 

a 25 

+ 

a 21 

a 22 

a 25 

= 2 

> 

a 21 

a 22 

a 25 


_ a 51 

a 52 

a 33_ 


_ a 51 

a 52 

a 33_ 



a -1 

L 3 1 

a 52 

a 33_ 


(4.4.6) 


(4.4.7) 


in which the 2 may be thought of as a coefficient of multiplicity, or, in other words, the coeffi^ 
cient which indicates how many equal matrices are to be added. If we adopt this notation, the 
following equation may be used to define it: 



a ll 

a 12 ••• 

a il^ 


a n 

a i2 ''' 

a lm 


ka n 

ka !2 

ka.. 

lm 

K 

a 21 

a 22 *•• 

a 2m 

- 

a 21 

a 22 

a 2m 

K = 

ka 21 

ka 22 

.. . ka 2m 


• 

• 

• 


• 

• 

• 


• 

• 

• 


a nl 

a n2 

a nm 


a nl 

u 

a n2 

a 

nm 


ka nl 

ka n2 

ka nm 


or in other words, the same constant may be "factored out" of every element of a matrix. 
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4.5 Multiplication of Matrices 

The product of an m x n matrix A and an n x p matrix 3 is the m x p matrix C defined 
by the equation 


'U 


E 

k=l 


a ik • b kJ ‘ 


(4.5.1) 


If the operations indicated by (4.5.1) are done by hand, it is convenient to organize the work 
as indicated in Fig. 4.5.1. In this figure an element of C is found by adding the products of 


a ll 

a 12 •" 

a ln 

a 21 

a 22 ''• 

a ln 

_ a ml 

a m2 "• 

a mn_ 


For example, the 

°12 

1—1 
r—l 

1* 

b 12 ••• 

b lp 

b 21 

b 22 •" 

b 2p 

_ b nl 

b n2 •" 

b np 

” c ll 

C 12 *•* 

C lp“ 

C 21 

c 22 •" 

°2P 

_ C ml. 

°m2 

c 

mp_ 


Fig. 4.5.1 


of intersection between row 1 of A and column 2 of 3 so its value is obtained as shown below 


J 12 


22 


Fig. 4.5.2 


i b n2 


a ll a 12 *‘* a ln 


C 12 a ll b 12 + a 12 b 22 + 


+ a ln b n2 
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Again, we may see from Pig. 4.5.1 that the process of generating each element of the product 
is only possible if the matrix used as the multiplier has as many columns as the matrix used 
as the multiplicand has rows. Therefore, the product of matrices is more restricted than the 
product of ordinary numbers. However, it should be obvious that the product of two square 
matrices of the same order is always possible. 


Example 4.5 

Find A-B if A 


■[: 3 • *■[; 


-2 

5 


p 

-2 

*1 

b 

5 

oj 

[2 5 ] R 2 + 15 ) 

(-4+25) 

(10+0 jl 

L 3j [(1+9) 

(-2+15) 

(5+0) J 


17 

10 


21 


15 


10 


The most common type of matrix multiplication encountered in these notes is the product 
of a square matrix and a column matrix or vector. In this case matrix multiplication is very 
simple since it becomes the sum of the products of the elements of each row of the matrix 
times the elements of the vector, as shown below for a 5 x 5 array. 


l ll a i2 a i3 


b ll 


a ll b ll + a 12 b 21 + a 13 b 31 

l 21 a 22 a 23 

• 

b 21 

= 

a 21 b ll + a 22 b 21 + a 23 b 31 

31 a 32 a 33_ 


_ b 31. 


_ a 31 b ll + a 32 b 21 + a 33 b 31_ 


Pig. 4.5.3 

The multiplication is simplified if the diagram shown in Pig. 4.5.4 is used. Where the brackets 
enclosing the vector indicate that the vector is really a column vector which has been written 
in row form. 


a ll 

a 12 

a 13 

a 21 

a 22 

a 23 

a 3i 

S 32 

a 33. 


b ll b 21 


b 3l}‘ 


= c 


Fig. 4.5.4 
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If Figure 4.5.4 is then written as shown in Fig. 4.5.5, it is easy to see that matrix multiplication 


becomes the sum of the product of each element of the matrix with its corresponding element 
the vector 


" a ll 

a 12 

a 13 


a n b n 

+ 

a 12 b 21 

+ 

a 15 b 31 

a 21 

a 22 

a 25 
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a 21 b ll 

+ 

a 22 b 21 

+ 

a 25 b 31 

_ a 51 

a 32 

a 53 



+ 

a 52 b 21 

+ 

a 33 b 31_ 


j b ll b 21 b 3l] 

Pig. 4 . 5.5 


4.6 Matrix Equations and the Ideograph 
Consider the equation 


V + a l it + a 2 = °- 

dt 


(4.6.1) 


We may define the operators 


3 = 3t 


s 2 = d 


dt 


, ... i 


dt 


(4.6.2) 


so that 


= sy * ^ = s 2 y* 

dt 


(4.6.5) 


and the equation may be rewritten in operator notation 


1 2 

a O y + a 1 s y + a 2 s y = 0. 


(4.6.4) 


If we make the coefficient of the highest derivative equal to 1, we obtain 


a 

a 


o 

2 


y 


1 2 

— sy + s y 
a 2 


0. 


(4.6.5) 


We will now show that this equation may be written in the matrix notation of (4.6.6), and that 
this notation will permit us to find an analog circuit which will yield the solution of the 
equation. 
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In this equation the matrix 


s -1 0 

O s -1 



is called the matrix of coefficients or coefficient matrix, the column matrix 


y 

sy 

_s 2 y_ 

is called the solution matrix or solution vector , and the matrix to the right of the equal sign 


0 

0 

0 


is called the matrix of constant terms . 

This equation yields, according to our definition of matrix multiplication. 


sy - sy + 0 


0 + 




s 2 y 


+ s 2 y 


"O' 

0 

0 


(4.6.7) 


The column matrices 


0 

a 

o 



at each side of the equal sign in (4.6.7) 

sy - sy + 0 = 0 

2 2 

+ sy-sy = 0 

a i 2 

y+ — sy+sy = 0 

a 2 


are equal if and only if 


Tnese equations do not contain any contradictions and, in fact, generate the given differen¬ 
tial equation. We will find that several matrix equations similar' to (4.6.6) may be found which 
have the property of generating the given differential equation. For example, the reader 
should verify that (4.6.9), (4.6.10) and (4.6.11) are three such systems. 
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(4.6.9) 


(4.6.10) 


(4.6.11) 


The reason for going through all the trouble of writing the differential equation in matrix 
form is that it has been shown (see references 1, 2, 3* 4) that it is possible to draw the 
analog circuit corresponding to the equation in such a way that there is a perfect correspon¬ 
dence between the matrix coefficients and the amplifier connections. This is done by drawing 
the analog circuit shown in Pig. 4.6.1. This figure, which is called an Ideograph , may be 
drawn directly from the matrix equation. 



In this figure there are three triangles representing amplifiers. These amplifiers are of the 
double-output type described in section 3.4. A detailed view of one of these amplifiers is 
given in Fig. 4.6.2. 
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Note: 

There is a change 
in sign between 
the input line and 
the POSITIVE out¬ 
put line 


Pig. 4.6.2 

The main diagonal of the coefficient matrix is used to determine whether an amplifier 

is an integrator or an adder. If the diagonal element contains the s operator, the amplifier 

is an integrator, otherwise, it is an adder, a constant multiplier, or a sign-changer. The 

non-diagonal entries of the matrix correspond to the gains of the amplifiers. Since there 

are four non-zero off-diagonal elements, four circles with appropriate gains are shown. 

The solution vector is translated Into the label of the outputs at the bottom of the 

0 

ideograph. Since double-output amplifiers are used, both y and -y, sy and -sy, and s y and 
2 

-s y are available for recording. Alternately, a single output amplifier may be used for 
simplifier 1, since only the normal positive output of this amplifier was needed. 

Since most computer Installations do not have double-output amplifiers, these may be 
simulated by a regular amplifier and a sign changer as shown in Pig. 4. 6 . 5 . 


Pig. 4.6.3 


As will be seen later, more efficient Ideographs may be developed for single-output 
amplifiers, so that the reader should not be too concerned about this detail now. 

The ideograph of Pig. 4.6.1 may be called the gain Ideograph to differentiate it from 
those developed in reference l 1 . It may be used very conveniently for making connections since 
it indicates very clearly the relationships between amplifier Inputs, outputs, and gains. For 
example, if we consider the output of each amplifier in the order of amplifier 



1 The gain ideograph was developed by the author as a natural extension of those developed 

by Honnell and Horn. It has the advantage of being less circuit oriented since it stresses 
the logical components of the analog circuit rather than its circuit components. 
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number, we find that the output of amplifier 1 is connected to one of the inputs 

of amplifier 5 with a gain of a^a^. The negative output of amplifier 2 is connected to the input 
of amplifier 1, with a gain of 1. The positive output is connected to the input of amplifier 
3 with a gain of a^/a^. The outputs of amplifier 3 are connected in the same manner. When all 
the outputs have been considered, all the inputs will be properly connected. Of course, the 
same result will be produced if the inputs, rather than the outputs, are connected in sequence. 

Some readers might prefer to draw the familiar block diagram corresponding to the ideograph 
shown in Fig. 4.6.1. This has been done in Fig. 4.6.4, It is somewhat inconvenient because 
we usually draw the block diagrams so that the derivatives decrease in power from left to 
right. The next section will show an ideograph which removes this inconvenience. 

Fig. 4.6.4 

We should also notice that an array which approximates the ideograph even more closely may be 
obtained by placing the solution vector below the matrix of coefficients els shown in Fig. 4.6.5. 




Fig. 4.6.5 


In this cstse the solution vector is enclosed in | j to warn the reader that such a transfor¬ 
mation has been applied. 
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Pig. 4.7.1 

When equation (4.7.1) is written in the form of (4.7.2) from which we may obtain Pig. 4.7.2, 
it is easy to see the correlation between the matrix equation and the ideograph of Pig. 4.7. 3 . 
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A question which might come to the reader’s mind is: What happens if the systems of (4.6.9), 
(4.6.10), or (4.6.11) are used as models for the ideograph? The answer is that each equival¬ 
ent matrix equation would yield an ideograph with perhaps a different appearance and, in general, 
a different number of amplifiers. 

The techniques developed in these notes will assume the use of a matrix of a particular 
form as the starting point of the procedure of drawing the ideograph. It will then be shown 
that one system may be converted to another by suitable matrix manipulations. 


4.8 Standard Form of Ordinary Differential Equations with Constant Coefficients 
In order to convert the differential equation 


ax + a n sx + a 0 s 2 x + ... + a^ -.s 11 ” 1 + s n = 0 
o 1 2 n-1 

into a matrix equation, we may use the standard form of (4.8.2). 


(4.8.1) 
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(4.8.2) 


In this system, all elements along the diagonal except the one in the lower right hand 
corner are equal to s. To the right of each s there is a -1, The last row consists of the 
coefficients of the differential equations. The column matrix of variables consists of all 
the derivatives of x in ascending order from top to bottom. 

Equation (4.8.2) leads to an ideograph in which the derivatives increase in order from 
left to right. An ideograph in which the derivatives decrease from left to right may be found 
by reflecting (4.8.2) to obtain(4.8.3). 
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In this system the operators appear in the other diagonal of the matrix, except the lower 
left hand position which contains a 1. To the left of each operator there is a -1, and the 
last row contains the coefficients of the differential equation. The solution vector consists 
of all the derivatives of x in descending order, from top to bottom. 

In both (4.8.2) and (4.8.3) the differential equation is in standard operator notation, 
that is, the coefficient of the highest derivative has been adjusted to unity. 

Most of the methods in this chapter are based on using (4.8.2) or (4.8.3) as the starting 
point and on transforming these equations into other equivalent systems by means of elementary 
matrix transformations. The advantage of these methods is that it is easier to manipulate 
matrices than to work with block diagrams, especially if this manipulation is being done by a 
digital computer program. 


4.9 Elementary Matrix Transformations and Minimization of Sign-Changers . 

Going back to matrix equations (4.6.6) and (4.7.1) we notice that a sign changer is needed 
for each negative entry in the coefficient matrix. This suggests the possibility of minimizing 
the number of sign changers by minimizing the number of negative matrix elements. Since, as 
was shown in section 4.6, there are many different matrix equations for each differential 
equation, it should be possible to find an equivalent system with fewer negative entries. One 
such system is given in (4. 9 . 1 ) which has been written with (4.6.6) for comparison. 

We can see that (4.9.1) may be obtained from (4.6.6) by changing the sign of all the 
elements in row 2 of the matrix equation and of all the elements in column 2 of the coefficient 
matrix. As a result of a double sign change, the diagonal element is not affected. 
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An operation such as this which does not affect the ability of the matrix equation to 
generate the given differential equation is called an "elementary matrix transformation." 

Elementary matrix transformations (EMT) are important because they permit us to transform 
one matrix equation to another. Since each equation represents a different analog system, we 
have at our disposal a method by which we may choose a system, from all those possible, which 
is superior in some respect to all others. For example, if we wish to minimize the number of 
double-output amplifiers, we should change the sign of those rows and columns which contribute . 
the largest number of negative entries in the coefficient matrix. Therefore, we may state the 
first elementary matrix transformation (EMT) as follows: 

EMT 1. - "The signs of all the elements in any row of the matrix equation may be changed if at 
the same time the signs of all elements of the corresponding column of the coefficient matrix 
are changed." 

Example 4.9.1 

The systems below are equivalent (EMT 1 on row 3 ) 



Example 4.9.2 

Minimize the number of sign-changers in the equation 



Solution: Applying EMT 1 to row-column 2 we obtain 
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Applying EMT 1 to row-column 4 we obtain 
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1. 


—1 

>5 
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03 

1 1 

3 of warning 

is 

now in 

order. 


0 


to a matrix in which the s operators appear along the main diagonal than to a reflected matrix. 
For example, if the matrix of Example 4.9.2 is reflected, we obtain 


1 

4- 1 y \ 

~ 0 0-1 s~ 

4- 

r s3y i 

z 

0 -1 s 0 


s 2 y 

3 

-1 s 0 0 

z. 

sy 

4. 

__1 -2 5 3_ 

1 

y ^ 


0 


(4.9.3) 


Notice that columns should be counted from right to left and rows from top to bottom in the coeffi¬ 
cient matrix. Rows in the solution vector should be counted from the bottom to the top. Changing 
the sign of row 3 may lead to error since the variable corresponding to row 3 does not line-up 
with row 5 of the coefficient matrix. 

The best way to apply EMT 1 to a reflected matrix is to place the solution vector below 
the coefficient matrix and to change the sign of the variable corresponding to the column 
which was changed in sign. 

In the matrix above we whould proceed as follows: 

Changing signs of row-column 3 gives. 
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Changing signs in row-column 1, 
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If this system is reflected, we obtain 
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which differs from the final result of Example 4.9.2 in that the signs of the elements of the 
solution vector are changed. This leads to EMT 2. 

EMT 2 

"Changing the sign of every element in the solution vector is an elementary matrix 
transformation." 

Other elementary matrix transformations which are useful in finding equivalent analog 
networks will be introduced as they are needed in succeeding sections. 

In order to show the use of (4.8.2), the following example is given. 


Example 4.9.5 Solve the equation 



The matrix equation is then constructed: 

(-? .| S -|. 2 + » 3 )-y.o. 

The matrix equation is then constructed: 



The use of EMT 1 on row-column 3 gives 



(4.9.5) 


(4.9.6) 


(4.9.7) 


which should be verified by the reader as an equivalent system by performing the matrix 
multiplication. 

Using EMT 1 again on row-column 1 we obtain 
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The corresponding ideograph is given in Fig. 4.9.1 
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Fig. 4.9.1 



Fig. 4.9.2 


The reflected gain ideograph of Fig. 


4.9.2 may be obtained by using the mirrow image of (4.9.8) 


0 0 1 s ■ 
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0 1 3 0 


-s 2 y 

1 s 0 0 
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or the mirror image of the ideograph of Fig. 


= o 


4 . 9.1 


( 4 . 9 . 9 ) 


4.10 Strategies for the Reduction of Negative Entries 

Whenever the number of variables is large and the coefficient matrix has many negative 
entries, it is not simple to decide what to do in order to minimize the number of negative 
entries unless some strategy is followed. We will show how we may do this by considering 
a larger matrix equation. 
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(4.10.1) 
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The first 

step 

is to prepare 

a "sign map" in 

which only the sign of the non-zero, non- 

diagonal entries 

, are 

shown. 

Also 

, the 

row numbers are carried on the right side of the map 
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0 -3 -3 -1 -3 -2 

pig. 4.10.1 

The second step is to add the number of positive signs and subtract the number of negative 
signs of each row-column pair and to write these differences at the bottom of each column. For 
example, row-column 1 has a difference of 0, since there is one negative entry and one positive 
entry. On the other hand, row-column 4 has one positive and two negative for a difference of -1. 

The third step is to change the sign of the row-column pair which contains the largest 
negative difference. Sometimes, as is the case in this example, there is more than one poss¬ 
ibility. In this case any of the choices will produce the same final result. 

Changing the sign of the second row and column we obtain 


0^ + 
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0 0 

0 0 

0 0 

+ + 
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i 


2 5 -1 -1 -3 0 

Fig. 4.10.2 

Notice that the sign of the row number was changed in order to remind us that the sign of the 
corresponding variable should be changed in the transformed matrix equation. 

The new sign map has a maximum negative difference In column 5, so we proceed to change 
the sign of the fifth row and column. If we proceed this way, the two following consecutive maps 
are produced. 
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Pig. 4.10.5 

As soon as all the differences are positive It Is Impossible to reduce the number of 
negative entries any further. We then use the sign map in order to obtain the matrix equation 
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(4.10.2) 


which requires a single sign-changer. 

The strategy which we have used in the preceeding procedure is essentially an algorithm 
which may be programmed on a digital computer. The digital computer would accept the sign map 
of Fig. 4.10.1 and change the sign of columns and rows until the column differences are all 
positive. 


4.11 Incorporation of Sign-Changers into the Ideograph 

Suppose that a matrix equation has been transformed by a number of elementary transfor¬ 
mations and we find that one or more entries are still negative. For example, the matrix equation 


(4.11.1) 
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(4.11.2) 
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by suitable transformations. Any attempts to eliminate the remaining negative entry will 
bring no improvements. This means that at least one amplifier will need a sign-changer associ- 

p 

ated with it, namely the amplifier which yields -s x. This is shown in the gain ideograph of 
Fig. 4.11.1. 



Fig. 4.11.1 

A very simple procedure may be used for obtaining an ideograph in which all the amplifiers are 
of the single-output type and all the sign changers are incorporated in the ideograph. The 
procedure will be presented by applying it to equation (4.11.2) 

1) By the use of elementary matrix transformations or by using the sign map, reduce the 
number of negative entries to a minimum (this has already been done in this example). 

2) Augment the coefficient matrix by adding a new row-column pair. (in this example, row 
and column 5). The elements of this new pair are: 

a) 1 in the diagonal position (5*5) 

b) 1 in the column of the negative constant (5*3) 

c) the negative constant, with its sign changed, in the row 
of the negative constant (2,5) 

(If there is more than one negative constant In the 
same column , repeat step c) 

d) zeros elsewhere. 
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5 ) Place a zero In the position occupied by the negative constant(s). 

4) Augment the column matrix by adding a row corresponding to the new row in the coefficient 

matrix (row 5)* Place in this position the negative of the variable in the row of the column 
matrix corresponding to the column of negative constant(s) in the coefficient matrix. 

s 1 0 

0 s 

0 0s 
a Q a x a 2 

_o 0 1 

x -sx (j-s x 

Actually, what we have done by following these steps is to add a new node to the system, 
a node being represented by an amplifier. This shows the power of the matrix notation, and it 


explains why 

there 

are 

many possible computer connections for each differential equation. Since 

the matrix may be 

transformed 

by a 

set of elementary operations, many alternate solutions may 

be obtained. 

For example, if 

the 

P ”Z 

node is added to (4.11.2) between -s x and s^x, we obtain 
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(4.11.5) 
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The gain ideograph of (4.11.5) is given in Pig. 4.11.2 
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Pig. 4.11.2 

The example below will illustrate the procedure of this section when a reflected equation 
is used as the starting point. 


Example 4.11 

Given the differential equation 



find the reflected ideograph. 

Solution: The equation is written in standard operator notation 
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Both row-column pairs 1 and 2 have sign sums of -2. 
Applying EMT 1 to row-column 1 we obtain 


Row-columns 5 and 4 have sign sums of -1. 
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EMT 1 on row-column 3 gives 



js^p -s^p sp -p 
Adding a column and a row 



The ideograph shown in Fig. 4.11.3 follows Immediately. 



Fig. 4.11.3 
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4.12 Magnitude Scaling by Matrix Manipulations 

Magnitude scaling may be accomplished by applying elementary matrix transformations to 
the matrix equations. Two more transformations which are needed will first be developed. 
Notice that the equation 


(4.12.1) 
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may be shown to be true by matrix multiplication. This leads to EMT 3 . 


EMT 3 

"An element in the jth row of the solution vector may be multiplied by a constant if 
every element of the jth column of the coefficient matrix is divided by the same constant'.’ 
Also notice that an equation is not affected if both sides of the equation are multiplied 

by the same const suit. The equations 
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are also equivalent. This may be stated: 


, etc. 


EMT 4 


"A row of the coefficient matrix may be multiplied by a constant if the corresponding 
row of the matrix of constants is multiplied by the same constant." 
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Magnitude scaling may now be shown wl-th the aid of an example. Example 2.5.1, presented 
in section 2.5 will be used for this purpose. 

The equation was 

x(t) + 2x(t) + 5x(t) = 0 

x(0) = 20 in (4.12.2) 

x(o) = 20 in/sec 

As we recall, the estimated maximum values dictated that the display variables be —, x/2 and 
5x. This means that these desired solutions must appear in the solution vector of the matrix 
equation rather than the usual variables x, x, and x. This may be accomplished as follows: 
l) Write the standard matrix equation for the differential equation 


(4.12.5) 


1 12 

2) Introduce the factor of 5 into x, the factor of into sx, and the factor into s x by 
EMT 5. 


s 

-1 

0" 


X 

0 

s 

-1 

• 

sx 

J> 

2 



A 


| -2 0 - 


~5x 

0 2s -5 

• 

SX 

_1 4 3_ 


1 S 2 X 
s X 


(4.12.4) 


5) Since the gain ideograph requires that the coefficients of the diagonal elements be equal 
to unity, EMT 4 is used on each row to secure this. Row 1 is multiplied by 5* row 2 is divid¬ 
ed by 2, and row 5 is divided by J>. 


T 1 

1 

H* 

O 

0“ 


5x 


5*0 ~ 

0 s 

.1 

2 


*1 SX 

= 

0 

2 

1 4 

J 5 

1 


j S 2 X 


0 

5 


(4.12.5) 


This is the scaled equation. Application of EMT 1 to row-column 2 and addition of a node 
gives the system 


0 

0 

iH 

03 

1_ 


X 

0 s | 0 


- -^SX 

1 _ . 4 


1 2 

3 0 1 3 


5 s x 

_° ! ° 1 _ 


_~5 3X _ 


O 


(4.12.6) 


The ideograph shown in Fig. 4.12.1, which results, is in effect the reflected form of the one 
obtained by using the highest derivative method. 
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Pig. 4.12.1 

Another way to approach magnitude scaling is to write the differential equation in terms 
of the new display variables as we had done before: 


(40* ! -0 

(5x o ) = 100 

(^2l) = 10 

and to use the matrix equation 


s -K 1 0 


’ (5x)“ 

0 s -K 2 

« 

■| sx 

1 4 1 


1 n 2 v 

J ? 


S X 


(4.12.7) 


(4.12.8) 


where K 1 and K 2 are adjusted to satisfy the first two rows of the product matrix. 
This requires that 


5 s x - K x - =0 


1 

2 


s 2 x 



(4.12.9) 


\ K 1 = 10 

) K 2 =! 
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It is easy to see that in this case the -1 terms in the standard matrix equation now become 
constants which depend on the manner in which the display variables are defined. In general, 
if -Kj is such a constant in the Jth row and Vj is the corresponding display variable, then 

K, = , (4.12.13) 

J j+n 


where n is the number of simultaneous differential equations represented by the matrix equation 
(in this case n=l). Time scaling may be handled as before, by multiplying the gains of the inte 
grators by an appropriate constant factor. The integrator gains are given by (4.12.13), so that 
in general, this equation may be written 


K = t_ !Li 

J *r V J+ n 

where ^ is computer time and t problem time. 


(4.12.14) 


4.13 Node Elimination 

A technique which is very useful in network analysis and synthesis is the modification of 
a network by eliminating nodes. We will introduce the technique by showing how a node may be 
eliminated from the ideograph by operating on the matrix equation. 

A node may be eliminated from an ideograph by applying the following rule: 

"A node represented by the diagonal element a^ may be eliminated by removing the kth 
column and the kth row from the matrix equation and by replacing each remaining matrix element 
a^ in the coefficient matrix with the expression a^ _ a lk a kj M . 


Suppose that in solving the equation 
a^X + ajSX - a Q X = 0 


(4.13.1) 


2 

we do not wish to obtain S X in the analog network. 

If the equation is written in standard operator notation, the matrix equation of (4.13.2) 
results. 


r s 

-1 

0 


” x ~ 

0 

s 

-1 


sx 

a 

0 

a l 

1 

« 

A 

OJ 

cd 

a 2 





(4.13.2) 


2 

Since we only wish to display x and sx, the node corresponding to s x may be eliminated 
by deleting row 3 and column 3 and by modifying each remaining matrix element according to the 
node elimination rule. This gives : 


83 


Digitized by 


Gck igle 


Original from 

UNIVERSITY OF MICHIGAN 



Generated on 2015-10-14 00:36 GMT / http://hdl.handle.net/2027/mdp.39015040284898 
Public Domain, Google-digitized / http://www.hathitrust. 0 rg/access_use#pd-g 00 gle 




-sx + sx = 0 

- t 2 * + (3 + -i) sx = - -2 X + s 2 x + -i sx = 0 (4.13.6) 

“2 a 2 a 2 a 2 

which reproduces the given differential equation. Elimination of node 5 yields an unfamiliar 
^1 

term (s + —) as one of the diagonal elements. It can be shown that this term may be interpreted 
2 a l 

as adding an extra gain equal to — across amplifier 1 as shown in the ideograph of Pig. 4.15.1. 

Pig. 4.15.1 

This indicates that the output of integrator No. 2 must be connected to its input through a 
gain of a 1 /a 2 

We conclude that elimination of node 5 has produced a simpler circuit which is preferable 
whenever the second derivative is not needed. 

Node elimination, however, has its practical limitations when applied to analog networks. 
This is because diagonal terms of the form (s + k) may be interpreted as was done above, but 
terms of the form (s - k) may not. Although in some cases it is possible to change the sign 
of the constant term by elementary matrix transformations, in some cases this may not be 
possible. For example, in the system. 

84 



Digitized fr. 


Gougle 


Original from 

UNIVERSITY C£|1)CHJ(aW^ 



Generated on 2015-10-14 00:36 GMT / http://hdl.handle.net/2027/mdp.39015040284898 
Public Domain, Google-digitized / http://www.hathitrust. 0 rg/access_use#pd-g 00 gle 


0 


(n.13.7) 




- 



s 

-1 

0 


X 

0 

s 

-1 

• 

sx 

a 

-a-, 

1 


S 2 X 

0 

1 

_ 




Elimination of node 3 gives 


“ 

- 



s 

-1 

# 

X 


s-a. 


_sx 

L 0 

1 J 




= 0 


The two negative constants may be eliminated by adding a node 


1 

a. 


x 

sx 

-sx 


(11.13.8) 


(11.13.9) 


but the procedure gives an ideograph with as many amplifiers as would be obtained if the second 
derivative was retained. 

This example has shown two things which we had not encountered before. One is the fact 
that when a node is added, all of the negative constants in the same column may be eliminated 
by moving them to the new column. Secondly, the negative portion of a term of the form (s - k) 
may be treated as any single negative constant and may also be moved, by itself, to the new 
column. 


4.14 Other Elementary Matrix Transformations 

Two additional elementary matrix transformations which are very useful in simplifying 
analog networks will now be given. The first transformation is stated as follows: 


EMT 5 


"A row of the coefficient matrix may be multiplied by a constant and added to any 
other row" 

Example 11.14.1 

The equation 


2 

-ax + a n sx + a 0 s x = 0 
o 1 2 

is satisfied by the matrix equation 


s 

-1 

0“ 


X 

0 

s 

-1 


sx 

«°l 

_L 


1 


s 2 x 

CVJ 

a 

_1 

a 2 





(4.14.1) 


(4.14.2) 
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as well as by the equation 

O (4.14.5) 


which results when row 5 is multiplied by 1 and added to row 2. This may be verified by perform¬ 
ing the matrix multiplication. 

Another elementary matrix transformation may be discovered if EMT 1 is applied on row- 
column 1 of (4.14.5) to obtain 


= 0 (4.14.4) 


If (4.14.4) Is now compared to (4.15.5), we notice that node 5 is actually not needed, since 
(4.15.5) satisfies the given differential equation and requires one less amplifier. 

Another way to reach the same conclusion is to notice that if node 5 is eliminated from 
(4.14.4) the remaining entries will not be affected, since all the elements of column 5 are 
zero except the diagonal. This leads to the concept of a redundant node , and to the following EMT. 
EMT 6 

’’Whenever all of the off-diagonal entries in a column are zero, the corresponding node 
is redundant and may be eliminated without affecting the rest of the matrix." 

EMT 6 suggests a method for eliminating nodes. The method consists on using EMT 5 to zero out 
all the off-diagonal column entries corresponding to a node, and then to apply EMT 6 to eliminate 
the node. 




4.15 Ordinary Differential Equations with Time-Dependent Forcing Functions . 
Consider the differential equation 

(a o s° + a 1 s 1 + a^s 2 + ... + a n s n ) x = f(t). 


(4.15.1) 


where the right side of the equation is a function of the independent variable t. 
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If the equation is written in the form 



We find that the matrix equation 


s 

-1 

0 . 

. . 0 

0 

0" 


X 

0 

s 

-1 . 

. . 0 

0 

0 


sx 

0 

0 

s . 

. . 0 

0 

0 


2 

S X 

• 

• 

• 

• 



• 

• 

0 

0 

0 

s 

-1 

0 


n-1 

S X 








n 

a 

a. 

a 0 

a ~i 

1 

-1 


s x 

0 

_1 

2 

n-1 




a 

a 

a 

a 





n 

n 

n 

n 





^0 

0 

0 

0 

0 

0 _ 


s f ( fc ) 







l n -* 


(4.15.2) 


(4.15.5) 


will satisfy the given differential equation. This matrix equation is somewhat unusual because 
all the elements in one of the rows are zero, including the diagonal. We will find, however, 
that this matrix will behave properly when applied any of the elementary matrix tranformations 
which we have used in connection with other matrices. 

Notice that node 6 is not redundant, since its diagonal element is not equal to 1. 

Example 4.15 

The equation 

a o y + a l + a 2 = f (t) (4.15.4) 

dt 

may be written In the form 


~s 

0 



Lo 0 0 oj 

{y sy s 2 y 4^} 


(4.15.5) 
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EMT 1 applied to row-columns 2 and 4 yields 


0 

1 


—— - TT— 1 


{’ 


-sy s y 


0. 

_ f(t) 


(4.15.6) 


The negative entry may be eliminated by adding a node to obtain 


s 1 

0 s 

^ 0 


I* 


0 

1 

-sy 


0 
0 

s^y 


0 

0 


2 V . Hi) 


0 “ 
0 

fl 

a 2 

0 

1_ 

sy 


(4.15.7) 


This ideograph is given in Pig. 4.15.1 



in^ot ta 

fw/icf • 
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Notice that the zero diagonal in row 4 indicates that node four does not contain an amplifier. 


Also notice that since the term -f(t)/a 2 appears as an entry in the solution vector, it will 
appear sis an output and will be represented as a vertical line in the ideograph. The actual 
function -f(t)/a 2 is the output of a function generator which may be placed in the location 
usually occupied by the amplifier. It is convenient to think of the horizontal line representing 
row 4 as the input line of the function generator. For example, if the forcing function in 
equation (4.15.4) is a function of y and s y rather than t, the connection shown in Fig. 4.15.2 
between the output of amplifiers 1 and 5 and the function generator input line would indicate 
the arguments of the function. 



LZ 





> 








3 


0 



z 


0 






3 





y 

-s 

/ V 

7 -Vfy. 

Cl 

s y - 

7 


Fig. 4.15.2 


4 .16 Equation with Variable Coefficients 
In the differential equation 


y + 



(4.16.1) 


the coefficient of the first derivative ^ is a function of the dependent variable y. If 
we write the equation in the form 


(y 




(4.16.2) 
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and define y ^ = f(y, sy) then the equation may be assumed to be an ordinary differential 
equation with a forcing function which depends on y and The matrix equation has the form of 

(4.15.5) and is given by 


~s 

-1 

0 

0“ 



" y ’ 

0 

3 

-1 

0 



sy 

1 

-1 

1 

-1 



s 2 y 

_0 

0 

0 

0_ 



_ysy_ 

the 

transformed 

matrix equation 

s 

1 

0 

0 " 


- 

y 

0 

s 

1 

0 



-sy 

1 

1 

1 

1 



s 2 y 

_0 

0 

0 

0 > 



-ysy 


= 0 


(4.16.3) 


(4.16.4) 


and the ideograph shown in Pig. 4.16.1 



Pig. 4.16.1 


4.17 Simultaneous Differential Equations 

The techniques of the previous paragraphs may be extended to the solution of simultaneous 
differential equations. For example, let us consider the differential equations 


x(t) + 2a y(t) + b 2 x(t) 
y(t) - 2a x(t) + b 2 y(t) 


= 0 
= 0 . 


(4.17.1) 
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These equations may be written in standard operator notation 


2 2 2 
s x + 0* s y + 0* sx + 2a-sy + b x + 0-y 


0 


0-s 2 x + s 2 y - 2a-sx + 0-sy + Ox + b 2 -y = 0 


(4.17.2) 

in which each equation contains every derivative. Then we may write the equation 


"s 2 

0 


X 


0 

2a* s 


X 


OJ 

0 


m 

V 

0 

8 2 . 

* 

_y. 

-h 

_-2a* s 

0 


.y. 

+ 

0 

1 

OJ 

• 

,J ■ 

A 


(4.17.5) 


Factoring out the solution vector 


s 2 0 


CQ 

cd 

OJ 

O 

1_ 


b 2 0 

1 


X 


0 " 

0 a 2 _ 

+ 

2a- s 0 

+ 

OJ 

0 

_l 

. 1 

• 

1 ' 

_y _ 

= 

0 


(4.17.4) 


Factoring out the s operator from the first matrix and the s operator from the second we 
obtain 


s 2 + 


0 

-2a 


2a 

0 


s + 


(4.17.5) 


This is an equation in which the coefficients of the operators and the dependent variable turn 
out to be matrices rather than single numbers. After a little practice, the matrix equation 
may be obtained by inspection from the original equations. 

Equations (4.17.5) may then be written in standard matrix form 


”s 

0 

1 -1 

0 

1 0 

0" 

0 

s 

1 0 

-1 

1 0 

0 

- - 

— — 

r- 

-- 

. — — - 

_ - . 

0 

0 

• s 

1 

0 

, -1 

0 

0 

0 

1 0 

s 

[ 0 

-1 

>’ 

0 

r - - . 

1 0 

1 

2a 

1 1 

0 

.0 

OJ 

,0 

: -2a 

0 

! 0 

1_ 

I. 

y 

sx 

sy 

s 2 x 

2 

s y 


(4.17.6) 


Notice that this equation may be written directly from the given equations. The derivatives 
of x and y increase from left to right, as before, with derivatives of the same order following 
each other. 

That equation (4.17.6) does not contain any contradictions and reproduces the two given 
differential equations may be shown by preforming the indicated multiplication. The dashed 
lines are used to emphasize the fact that this matrix form is very similar to the one we used 
for a single differential equation, except in that the matrices 


0 

CQ 


\ 1 
i— 1 

0 , 


0 

, 0 


t 


, and 


0 s 


0 

1 

M 

t 


L° oj 
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are used in place of s, -1, and 0, and that the original coefficients a Q , a^, ... , a n -i* 1 are 
now the matrix coefficients of (4.17.5). 

As in the case of a single equation, (4.17.6) is one of the many forms in which the 
simultaneous equations may be written. All of the elementary matrix transformations which applied 
to the matrix equations in the previous sections also apply to (4.17.6). 

If a'feign map" is prepared for this system, it will be found that a minimal amplifier 
system is given by 


s 

0 

1 

0 

0 

cf 


X 

0 

s 

0 

1 

0 

0 


y 

0 

0 

s 

0 

1 

0 


-sx 

0 

0 

0 

s 

0 

1 


-sy 

b 2 

0 

0 

-2a 

1 

0 


s 2 x 

0 

b 2 

2a 

0 

0 

1_ 


1 

GQ 

J _ 


O 


(4.17.7) 


Addition of 

a node 

produces 

the ideograph 

equation 


“s 
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1 
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0 

0 

0 " 
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s 

0 

1 

0 

0 

0 


0 

0 
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0 
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0 

0 


0 

0 

0 

s 

0 

1 

0 


b 2 

0 

0 

0 

1 

0 

2a 


0 

b 2 

2a 

0 

0 

1 

0 


_0 

0 

0 

1 

0 

0 

1 

i- 

y 

-sx 

-sy 

2 

s X 

s 2 y 

syj 


(4.17.8) 


Another example will further illustrate the method. 
Example 4.17.1 

Solve the differential equations 

]y(t) + 254 y(t) - 1296 9 =0 
j 8o(t) - 54 y(t) + 5616 0=0 
with initial conditions 

) y(0) = 100 y(0) = 0 

) 9(0) = 1.47 9(0) = 0 

and the display variables 

( 2^0 ) ’ (~ 20 ) ' y 

( | ) , (29) , (509) 
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The equation ia written in standard operator notation 

s 2 y + 0*s 2 9 + 0*sy + 0*s9 + 2^4y - 12969 = 0 

2 2 

0*s y + a 9 + 0*sy + 0*s9 - 6.75y + 7029 = 0 


(4. 17 .II) 


and in matrix notation 


"l 

o" 

s 2 + 

‘o 

o] 

. s + 

234 

-1296 

• 

y" 

0 

1 _ 


_0 

°J 

VJI 

702 

( 

9 _ 


The standard matrix equation is 


= O (4.17.12) 
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s9 
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0 
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0 


s 2 y 

^ 6.75 
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0 

0 

0 

1_ 


1 

00 

ro 

CD 

L- 


= C (4.17.13) 


The display variables are Introduced Into the solution vector by application of EMT 3 
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0 

-20 
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s 

50 

0 

1 
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0 

20s 

0 

-250 
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0 
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0 
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2s 9 

254 

-25.92 

0 

0 

250 

0 


fso 

^6.75 

14.04 

0 

0 

0 

5 


> 

l- 


rQ (4.17.14) 


Now each row of the coefficient matrix may be multiplied or divided by a constant by using EMT 4 
so that all diagonal elements become unity. This gives. 
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sy 
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0 

-10 


2 . s9 

.936 
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1 
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sfjr 

250 

: 1.35 

2.808 

0 

0 

0 
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s 2 9 

-~5~ J 


= O (4.17.15) 
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A sign map would indicate that row-columns 2, 2, and 6 should be changed in sign. This yields the 
ideograph equation 


(4.17.16) 


j y -5oe -fg ^ u -¥) 

Notice that every row containing s in (4.17.16) has a gain which is of the order of ten times 
larger than unity. This is to be expected since the estimated undamped frequencies for the 
equations may be found to be c*> *-15 rad/sec and °° q 0,6 30 rad/sec. This suggests that all of 
the integrator gains be divided by 10 to slow down the solution to a frequency range at which 
the analog components respond more accurately. 

An alternative procedure for solving this example is to write the given differential 
equations in scaled form 


(-gjy + • 9?6(y) - .104(506) = 0 
( 5)" 1 - 55 (y)+ 2.81 (500) = 0 

j y(°) = 100 (li 21 ) = 


\ y( 0) = 100 = 0 

j 500(0) = 73.5 20(0) = 0. 

Equation 4.17.15 then follows directly If equation (4.12.15), with n = 2, Is used In obtaining 
the K constants. This gives 


3 V c 


k 4 = v: 


s(2 s9) 
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4.1 Problems 


In the expression A* B we may say that 3 is pre - multiplied by A and that A is post - 
multlplied by 3 . Show that in general a different result is obtained if A is pre-multiplied 

by ft than if A is post-multiplied by 3 . 

4.2 If A*B exists, does B*A exist? 

4.3 Find the product P*Q of the matrices below 



"a 

-1 

O' 


a 

p = 

0 

a 

-1 

Q = 

a 2 



-2 

1_ 


.a 3 . 


4.4 Find R • S and S • R if R and S are given by 



p 

0 

0 " 


~ a n 

a 12 

a 13~ 

R = 

0 

1 

0 

S = 

a 21 

a 22 

a 23 


_0 

0 

K_ 


. a 31 

a 32 

a 33_ 

Reoeat problem 

4.4 

with 





’l 

0 

0” 




R = 

0 

K 

0 





0 

0 

1 





and write down a rule in regard to pre-multiplying or post-multiplying a matrix by another 
matrix containing unity in all diagonals except one. 

4.6 Show that A *C + ft • C - (A -*■ &) • C 

that is, that the distributive law of multiplication applies when the variables are matrices. 

4.7 Solve the following problems by the matrix approach: 

a) Problem 2.1 

b) Problem 2.2 

c) Problem 2.6 

d) Problem 2.7 

e) Problem 2.8 
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4.8 Find the gain ideograph corresponding to the equations: 

j 400x(t) + 1500 y(t) - 300J(t) + 20x(t) = 0 

1 1200x(t) + 2500y(t) + 200x(t) - 100y(t) + 50 y(t) = O 

and the display variables: 

2x(t) , 5y(t) , ^ x(t) , y(t) , TG x(t) , ^ y(t) 

4.9 Obtain an ideograph for Problem 4.7a in which the third derivative x(t) is not displayed. 

4.10 Eliminate a(t) from the ideograph of Problem 4.7b. 

4.11 Find the integrator gains k^ k 2 , kand k^ in the scaled system below and recommend a 

change in time scale if it would seem required 


s 

0 

- K i 

0 

0 

0 


’lOOx ' 

0 

s 

0 

- K 2 

0 

0 


500y 

0 

0 

s 

0 

- K 3 

0 


lOsx 

0 

0 

0 

s 

0 

- k 4 

• 

50sy 

1.25 

-.26 

0 

-1 

1 

0 


S 2 X 

2 

.-2.5 

1.8 

1 

0 

0 

1 


2 

_s y 


4.12 Obtain the minimal amplifier circuit for the matrix equation of Problem 4.11 by using 
the sign map technique. Add enough nodes to produce an ideograph which uses single- 
output amplifiers. 
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Answer to Problems 


Exercise 1.1 


a) 

-15 

volts. 

b) +5v. 

c) .15 

d) 

5.5 e) 20v. f) -5E X g) 5 

h) 

.3 


i) 10t+30 

J) \ t 2 -5 

k) 

-y( t) 1) - 50 t-y(t) m) -y(t) 

n) 

.1 

0) 

initial = 5 

gain-. 5 

P) 

top gain = 1, bottom gain = A, initial = -Aw(o) 

q) 

1 

2 ' 

3 

•1, J 

2 
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2.5 a) Amplifier Output 

1 +100 

2 0 
3 -100 

b) 1 and 4 


c) Amplifier 1 would try to reach +200 volts and would become overloaded. 



( 222 

y x Z -sy -sx -sZ s y s x s Z -20 sx 
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